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ABSTRACT 

The hot, X-ray-emitting gas in superbubbles imparts energy and enriched material to the interstellar 
medium (ISM) and generates the hot ionized medium, the ISM's high-temperature component. The 
evolution of superbubble energy budgets is not well understood, however, and the processes respon- 
sible for enhanced X-ray emission in superbubbles remain a matter of debate. We present Chandra 
ACIS-S observations of two X-ray-bright superbubbles in the Large Magellanic Cloud (LMC), DEM 
L50 (N186) and DEM L152 (N44), with an emphasis on disentangling the true superbubble X-ray 
emission from non-related diffuse emission and determining the spatial origin and spectral variation 
of the X-ray emission. An examination of the superbubble energy budgets shows that on the order of 
50% of the X-ray emission comes from regions associated with supernova remnant (SNR) impacts. We 
find some evidence of mass-loading due to swept- up clouds and metallicity enrichment, but neither 
mechanism provides a significant contribution to the X-ray luminosities. We also find that one of the 
superbubbles, DEM L50, is likely not in collisional ionization equilibrium. We compare our observa- 
tions to the predictions of the standard Weaver et al. model and to 1-D hydrodynamic simulations 
including cavity supernova impacts on the shell walls. Our observations show that mass-loading due 
to thermal evaporation from the shell walls and SNR impacts are the dominant source of enhanced 
X-ray luminosities in superbubbles. These two processes should affect most superbubbles, and their 
contribution to the X-ray luminosity must be considered when determining the energy available for 
transport to the ISM. 

Subject headings: ISM: bubbles — X-rays: ISM — H II regions — supernova remnants — Magellanic 
Clouds — open clusters and associations: general 



1. INTRODUCTION 

The hot phase (T > 10 6 K) of the interstellar medium 
(ISM) is generated primarily by the mechanical power 
imparted by stellar winds from massive stars and super- 
novae (SNe). The clustered nature of massive star for- 
mation necessarily dictates that some stellar winds and 
core-collapse SNe will occur in low-density cavities exca- 
vated by other winds or SN e, thus creating super b ubbles . 

Theoretical work by iMac Low fc McCravl (|1988l ) 
demonstrated that superbubbles could be considered 
more powerful versions of stellar wind-blown bub- 
bles, w hich already had a well-developed theoretical 
mode l (|Dvsonl 119731: ICastor et all 119751: IWeaver et aT] 
119771). The most de tailed of these papers was that 
of IWeaver et al.1 (|1977l ). which has become the de-facto 
standard model explaining the structure and physics of 
wind-blown bubbles and superbubbles. 

Superbubbles play important roles in generating and 
maintaining the multi-phase ISM in star-forming galax- 
ies (see e.g. iNorman fc Ikcuchi 1989: lNorman fc Ferraral 
[1991 lOev fc Clarke! 11997ft . The influence of superbub- 
bles can extend beyond the thin gaseous disk, as some 
superbubbles vent their hot, metal-enriched pl asma into 
the g alactic halo, creating galactic fountains (jBregmanl 
1980) in normal spiral galaxies like the Milky Way 
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and feeding galactic superwinds in starburst galaxies 
(jHeckman et al.lll990ft . 

As the hot, collisionally ionized plasma in superbub- 
bles emits primarily by X-ray line emission and ther- 
mal bremsstrahlung, X-ray observations provide a crucial 
window on these energetic phenomena. However both 
X-ray and optical observational studies indicate that our 
current understanding of superbubble physics and ener- 
getics is incomplete. 

The earliest imaging X-ray observations of star- 
forming complexes in the Magellanic Clouds with the 
Einstein Observatory revealed that so me superbubbles 
were indeed detected as X-ra y sources (|Chu fc Mac Low! 
ll990HWang fc Helfandlll991f ). Surprisingly, the detected 
superbubbles typically had X-ray luminosities an or- 
der of magnitude greater than predicted from the stan- 
dard Weaver et al. model (we refer to these as X- 
ray-bright superbubbles). These X-ray-bright superbub- 
bles demonstrate that the processes responsible for X- 
ray emission in superbubbles are not yet fully under- 
stood. Not all superbubbles are over-luminous, how- 
ever. Other superbubbles in the Large Magellanic Cloud 
(LMC) are as faint or fainter than expe cted based on the 
Weaver et al. model (|Chu et al.lll995l : referred to as X- 
ray-dim superbubbles). Subsequent observations using 
the Rontqensatelli t (ROSAT) and ASCA X-ray Obser- 
vatories (see e.g. lCIiu_eialJ [19931; iMagnier et all 119961 : 
iDunne et al.ll200l"ft have confirmed these results. 

Solutions to the X-ray-bright superbubble problem fall 
into two categories; (1) replacing the standard Weaver 
et al. model at some or all epochs of superbubble 
growth, and (2) supplementing and extending the stan- 
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dard Weaver et al. model. 

A popular and elegant example of the first approach 
was pioneered bv lChu fc Mac Low! (|1990f) . who proposed 
that the impact of the blast waves from any SNe that 
occurred near the supcrbubble shell can shock-heat the 
shell to X-ray-emitting temperatures. Using 1-D mod- 
els employing the K ompaneets thin shell approximation 
(Kompaneets 1960) they demonstrated that this shock- 
heating leads to a large, but temporary, increase in the 
total X-ray luminosity from the superbubble. Evidence 
exists for supernova remnant (SNR) shocks within su- 
perbubbles, and SNR shell impacts are also the most 
plausible cause of the unusually hi gh shell expansion ve- 
locitie s seen in some superbubbles (|Chu et ai1ll994t lOevi 
1996b). Nevertheless, their precise contribution to the 
X-ray luminosity is uncertain. 

The most-promisi ng example of a so lution of the sec- 
ond kind is that of ISilich et all (|2001h . The standard 
Weaver et al. wind-blown bubble model takes no account 
of metal enrichment once SNe occur; the metal abun- 
dances in the hot plasma are assumed to be the same 
value at all times. For a collisionally ionized plasma in 
the temperature range T = 10 6 — 10 7 K the X-ray emis- 
sivity is roughly proportional to the metal abundance 
of the plasma, as line and rec ombination dominate over 
pure thermal bremsstrahlung. ISilich et all point out that 
if the plasma is enriched in abundance by a factor of 10 
or more by SN activity, then a similar enhancement in 
the soft X-ray luminosity should be expected. In this 
hypothesis the structure of the hot X-ray-emitting su- 
perbubble interior remains that given by the Weaver et 
al. model at all times, but its chemical enrichment must 
be taken into account self-consistently. This hypothesis 
remains to be tested - as yet there is no observational 
evidence for high metal enrichment in X-ray-bright su- 
perbubbles. 

Constraining the extent to which these, or any other, 
mechanisms contribute to the total X-ray emission of su- 
perbubbles requires higher quality X-ray data than was 
available with the earlier generations of X-ray observato- 
ries (Einstein, ROSAT, and ASCA). Potential discrimi- 
nants between the different models hinge on spatial vari- 
ations in the X-ray spectral properties of the plasma 
within superbubbles, which requires observations with 
an instrument combining high spatial resolution, moder- 
ate to high spectral resolution and high sensitivity. Both 
modern X-ray observatories, the Chandra X-ray Obser- 
vatory and the XMM-Newton, fulfill these requirements. 

We have obtained and analyzed Chandra observations 
of two X-ray-bright LMC superbubbles, with the Ad- 
vanced CCD Imaging Spectrometer (ACIS) instrument. 
T he target obje c ts are DEM L50 and DEM L152 from 
th dDavies et all (| 19761 ) Ha catalog of LMC Hn regions. 
DEM L152 is also part of the nebula cataloged as N44 by 
iHenizd (|1956H Other properties of these superbubbles 
are highly constrained with existing multi-wavelength 
observations, properties such as the host stellar popula- 
tion, shell expansion characteristics, local interstellar gas 
density and intrinsic absorption. We adopt a distance of 
50 kpc to the LMC (jFreedman et al.ll200lD . At this dis- 
tance an angular offset of 1" corresponds to a physical 
size of 0.24 pc. 

DEM L50 has not previously been the target of pointed 
X-ray observations. Optical and Hi observations are re- 



ported in lOev et all (|2002l ). It is a roughly elliptical bub- 
ble with major and minor axes of 9' and 7' respectively, 
associated with the superimposed SNR N186 D on its 
northern edge. It is unclear whether the SNR is phys- 
ically connected to the superbubble or whether this is 
a chance projection of separate objects along a line of 
sight. Otherwise DEM L50 is relatively isolated from 
other Hn nebulae and young stellar clusters. 

In contrast DEM L152 is part of the large and rela- 
tively crowded N44 complex, a set of shells and Hil re - 
gions covering ~ 25' of sky (jMeaburn &: Laspiaslll99ll ). 
Given the limited 8' field of view of individual Chandra 
ACIS CCD chips we choose to concentrate on the X-ray- 
bright features Shell 1, the South Bar, and Shell 3 that 
were identified in earlier ROSAT 'and A SCA observations 
(|Chu et allll993HMagnier et al.lll99(l . 

By exploiting Chandra's ~ 1" spatial resolution we 
separate the true superbubble X-ray emission from non- 
related diffuse X-ray emission in the immediate vicinity 
of each superbubble and investigate the degree of spatial 
variations in soft X-ray properties within these bubbles 
(see § |4j) . If the diffuse X-ray emission in these bubbles is 
very highly structured, arising in small-scale regions as- 
sociated with the optical shells, then these Chandra ob- 
servation have the spatial resolution to detect this. In §0 
we compare o ur results to the predictions of the standard 
iWeaver et all model and to 1-D hydrodynamical simula- 
tions including SN-shell impacts and discuss superbubble 
energy budgets in light of our observations. We discuss 
the metallicity-enhancement scenario and the possibility 
of mass-loading from clouds in light of our observations 
in § M 



2. OBSERVATIONS 

DEM L152 was observed with the Chandra ACIS-S in- 
strument on 2002 September 22 for ~ 20 ks (observation 
identification number [ObsID] 3356), and DEM L50 was 
observed for - 40 ks on 2003 January 01 (ObsID 3355, 
see Table. [T]). The observations were made to accommo- 
date the majority of the optical nebula of each bubble on 
the more sensitive back-illuminated ACIS S3 CCD chip, 
but parts of each nebula spill onto the adjacent S2 and 
S4 CCD chips. 

Data reduction and analysis were performed using the 
Chandra CIAO software package (version 3.4) with the 
associated calibration database (CALDB version 3.3). In 
addition HEASOFT (version 6.1.1) was used for some 
tasks, including spectral fitting with XSPEC (version 
12.5.1). 

Each raw data set was reprocessed to take account of 
the latest calibration of important effects such as time- 
dependent CCD gain, charge transfer inefficiency and 
contamination of the optical-blocking filter. Periods dur- 
ing the observation that experienced higher-than-normal 
levels of Solar X-ray and particle events (commonly re- 
ferred to as background flares) were rem oved using the it- 
erativ e 3<7 clipping method described in lStrickland etall 
(2001). The effective on-source exposure time remain- 
ing in each ACIS CCD chip after this procedure is given 
in Table [TJ Only a small fraction, < 1%, of the total 
exposure time was lost to background flares in either ob- 
servation. This is much lower than the typical ~ 20% 
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TABLE 1 

Chandra observations OF DEM L50 and DEM L152 



Target R.A. Decl. Date Instrument ObsID Exposure time (CCD chip) 

DEM L50 04 h 59 m 47 s -70°11'37" 2003-01-01 ACIS-S 3355 37011.9 (S3), 36720.9 (S2), 36714.5 (S4) 

DEM L152 05' i 22 m 17 s -67°56'38" 2002-09-22 ACIS-S 3356 18713.6 (S3), 18713.6 (S2), 18713.6 (S4) 



Note. — The right ascension and declination values quoted (J2000.0 coordinates) are the approximate center 
of the superbubble cavity as seen in Ha images. The ObsID is the unique observation identification number 
assigned by the Chandra Science Center. The exposure times (in seconds) quoted are those remaining after 
the removal of periods of higher-than-normal background. The specific ACIS CCD chip to which the exposure 
time applies is noted in parentheses. 



fraction of ACIS observation time affected by flares 0. 

We used the wavelet-based source detection algorithm 
wavdetect to search for point-like objects in each X-ray 
observation. Each chip was treated separately, searching 
for sources in images created in the soft X-ray 0.3 - 2.0 
keV energy-band, hard X-ray 2.0 - 8.0 keV energy band 
and total 0.3 - 8.0 keV energy band. Only sources with 
signal-to-noise ratios > 2 were accepted as point sources 
in the initial iteration of the data analysis. In a few cases, 
low S/N features that most probably are genuine point 
sources were missed by the source-detection algorithm. 
These sources became apparent only after all brighter 
point sources had been removed. We assessed each such 
feature, and removed those that appeared to be point- 
like based on our personal scientific judgment. The point 
sources are shown marked on soft (E = 0.3 - 2 keV) and 
hard (E = 2 - 8 keV) energy band images of DEM L50 
and DEM L152 in FiglTJ 

We screened out the events from detected point sources 
in the data used for both imaging and spectral analysis 
of the superbubbles to a radius equivalent to 4 Gaussian 
a, based on a fit to the radius of the Chandra PSF as 
a function of off-axis angle. The holes left in any image 
by source removal were filled in using the CIAO task 
dmfilth. The observed distribution of pixel values in a 
background annulus of thickness 6 ACIS pixels (~ 3") 
around the point source defines the probability density 
function from which random values were chosen to fill 
in the source region (The POISSON option in dmfilth). 
Care was taken to ensure that the background annulus 
chosen did not contain other point sources which would 
bias the interpolation. 

For spectral analysis, sources detected in any of the 
three bands (the soft, hard and total bands described 
above) were excluded when creating diffuse-emission- 
only spectra. The soft band source lists were used for 
any image of the diffuse emission created within the en- 
ergy range 0.3 - 2.8 keV. The hard and total band source 
lists were used for any image created within the 2.8-8.0 
keV energy ban d. We discuss the data reduction of the 
images in § 13.11 and the data reduction of the spectra in 
§1121 

3. IMAGE ANALYSIS 

3.1. Image Data Reduction 

For the purposes of background subtraction we used 
the blank-sky datasets provided as part of the Chan- 

5 See "General discussion of the quiescent and flare compo- 
nents of the ACIS background", by Maxim Markevitch (2001) at 
http : //exc .harvard. edu/contrib/maxim/bg/index .html 



dra CALDB, as the diffuse soft X-ray emission from 
both DEM L50 and DEM L152 covers a large fraction 
of the S3 chip in each observation. These blank-sky data 
sets comprise hundreds of kiloseconds of data per CCD 
chip, providing a high S/N estimate of the typical back- 
ground experienced by the ACIS detectors. Since the 
net background experienced in the ACIS detector does 
change with time, the blank sky data was re-normalized 
so that the mean X-ray surface brightness (excluding 
point sources) for each chip in the E = 3 - 6 keV en- 
ergy band matched that in each observation. Typical 
variations in the hard X-ray background in the ACIS in- 
strument are ~ 2% (root mean square) from observation 
to observation (jHickox fc Markevitchl [2006h - 

We used a hard X-ray band to calculate the renormal- 
ization factor for the blank-sky data sets as soft diffuse 
X-ray emission is present in both the S2 and S3 chips 
in the two observations (see Fig. [1]). The majority of 
the superbubble emission falls within the field of view 
of the S3 chip in both observations. Part of DEM L50 
spills ove r onto the S2 chi p, however, and N44 shell 3 (a 
SNR; see iChu et all 119931 ) is a strong soft X-ray source 
in the S2 chip of the DEM L152 observation. Excess soft 
X-ray emission is not apparent on the S4 chip in either 
observation. 

Note that apparently diffuse hard X-ray emission in 
a s uperbubble has only been detected in 30 D oradus 
C (ISmith fc Wand [2001 lYamaeuchi et all l20Tol) . Our 
method of re-normalizing the blank-sky data sets means 
that there will be no significant diffuse X-ray emission in 
the E = 3 - 6 keV energy band, once point sources and 
the background have been subtracted. This approach 
is valid, in that a visual inspection of the raw, non- 
background-subtracted data reveals no evidence for hard 
X-ray emission within DEM L50 or DEM 152, and the 
mean source-free hard X-ray surface brightness in each 
chip of our observations is less than or equal to that in 
the blank-sky data. We cannot rule out the presence of 
diffuse hard X-ray emission in our observations of DEM 
L50 and DEM L152, but if such emission is present, it is 
very faint. 

For DEM L50, the S3 chip needed essentially no renor- 
malization, as the observed source-free E = 3 — 6 keV 
surface brightness differed from the blank-sky data sets 
by only 0.9 ± 1.8%. For the front illuminated S2 and 
S4 chips the blank sky backgrounds were normalized 
downward by 5.9% (uncertainty 2.1%) and 14.9% (un- 
certainty 2.0%) respectively. It is not surprising that 
the hard background surface brightness in the blank-sky 
data might be higher than in our observations of DEM 



4 



Jaskot, Strickland, Oey, Chu & Garcfa-Segura 



10 - 



t 

o 

L. 

o 

LJ 
Q 
■■:] 



O - 



in k 



(a) 



«3 



>©' 01 

'o- 
a 



1 



Qo - 



in 



(b)\ 



~l 1 1 1 1- 



* • i 



r — 1 — r 

i 

i 



\ 

\ 

1 



19 ^--^"'l 



-= 

I 

_4 



■ ■ i 



-5 







Q 

.0 O - 



o 

Q 
<! 



m 3 



(c) 



>0 O l 

o 



o g \ 



S4 



S3 







S2 



i — i — i — i — i— p — . — r 
i 



' • • . i S2 



I 



'S3 

. I 

\ 

1 

I 

• • ■ 1 



o 



o 

i_Q_ 



-5 

ARA (arcmin) 



-5 

ARA (arcmin) 



Fig. 1. — X-ray emission in the vicinity of DEM L50 (a and c) and DEM L152 (b and d). Each image shows the raw, unsmoothed, photon 
distribution in the soft 0.3-2.0 keV (panels a and b) and and hard 2.0-8.0 keV energy bands (panels c and d), binned in O'.'984-wide pixels, 
over a ~ 240 X 240 pc (16'4 X 16f4) region centered on the coordinates given in Table[l] Point-like X-ray sources detected in the soft and 
hard energy bands are shown surrounded by a circle, equivalent in size to the region used to remove the point sources from the images and 
spectra. The images are displayed using a asinh intensity scale. 



L50 and DEM L152. The blank sky data sets are con- 
structed from a large number of data sets and thus from 
data experiencing a higher (but more normal) fraction of 
mild background flares than our observations. 

For the DEM LI 52 observation the blank sky back- 
ground was normalized down by 4.4% (uncertainty 
±2.4%). For the S2 and S4 chips the blank sky back- 
grounds were normalized downward by 2.0% (uncertainty 
3.0%) and 16.6% (uncertainty 2.8%) respectively. 

The magnitude of the difference in the mean source- 
free X-ray surface brightness in the E = 3 - 6 keV energy 
band for the S4 chip in both observations is somewhat 
surprising. This CCD chip is known to suffer from sig- 
nificant flaws in its serial readout that leads to artifi- 
cial events being registered (visible as streaks in images) 
that can only be partially screened-out in software pro- 
cessing. We speculate that this effect is normally ex- 
acerbated by the unwanted particle events common in 



background flares, thus accounting for the pronounced 
difference between the blank sky data for the S4 chip 
and our observations. 

We altered the coordinate system of the blank-sky data 
to match the point of the appropriate observation. Thus 
our background subtraction for image analysis consisted 
of subtracting an image created from the blank-sky data 
from an image from the real observational data, scaled 
by the ratio of total exposure times in the observation to 
the blank-sky data for the appropriate chip. 

In Figures [2] and [3] we present adaptively smoothed 
X-ray images of DEM L50 and DEM L152 in a va- 
riety of soft X-ray energy bands, along with optical 
Ha+[Nn] imaging and Hi column density maps. Fig. [4] 
presents three-color composite images of the super- 
bubbles that combine narrow-band optical Ha+[Nn], 
[Oin] and soft X-ray emission. 

Adaptive smoothing using the CIAO tool CSMOOTH 
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Fig. 2. — Adaptively-smoothed background-subtracted images of the diffuse X-ray emission associated with DEM L50. Panel (a) shows 
the E = 0.3 - 2.8 keV emission on a square root intensity scale, overlaid with contours increasing in factors of two beginning at a surface 
brightness of T,x = 4 X 10 — 7 counts s — 1 arcsec - 2 . These contours are shown overlaid on a narrow-band Ha image from Oey et al. (2002) 
in (b). Panel (c) shows the E = 0.3 - 2.8 keV emission in grey scale overlaid with contours of Hi column density from lOev et al.l (2002) 



linearly spaced between Nn = 0.5 X 10 cm 



and 2 X 10 21 cm -2 . Panels (d), (e) and (f) show adaptively smoothed X-ray images in the 
E = 0.3 - 0.6, E = 0.6 - 1.6 and E = 1.6 - 2.8 keV energy bands respectively, using a square root intensity scale. Contour levels increase 
in factors of two beginning at a surface brightness of Sx = 2 X 10 -7 counts s -1 arcsec -2 . Dashed lines denote the edges of the Chandra 
ACIS CCD chips, and the dashed red lines in (a) and (f) show the spectral extraction regions. 



(jEbeling et al.l I2006T I attempts to differentially smooth 
an input image to achieve a relatively uniform local S/N 
ratio, so that bright regions are smoothed with a smaller 
smoothing kernel than low surface brightness regions. 
We initially smoothed the diffuse image in the E = 0.3 
- 2.8 keV energy band, after point source removal but 
prior to background subtraction, with a target S/N ra- 
tio of 3 but with a maximum smoothing kernel equiv- 
alent to a Gaussian of FWHM= 40". The smoothing 
map calculated in this step was then applied to all the 
background-subtracted soft X-ray diffuse images, so that 
any differences between images in different energy bands 
are intrinsic and are not due to the application of differ- 
ent smoothing maps. 

Given claims that CSMOOTH can generate spu- 
rious structure in adaptively smoothed images 
(|Diehl fc Statler! I2006D . we also investigated adaptively 
smoothed images using the completely independent 
adaptive smoothing algorithm implemented in the 
XMM-Newton data analysis software SAS as ASMOOTH. 
This algorithm generated images very similar to those 
produced by CSMOOTH, except that the stated statistical 
significance of individual features within the images 
differed significantly between the two methods. We 
also generated images using uniform smoothing with 
a 2-dimensional Gaussian mask of FWHM= 30" (not 
shown) in order to conservatively verify the existence of 
features seen in the adaptively smoothed images. 

3.2. DEM L50 



The diffuse soft X-ray emission in DEM L50 is clearly 
confined within the optical superbubble and SNR shells 
(Fig. [5b). The apparent X-ray surface brightness of the 
SNR N186 D (jOev et al J 120021 ) on the northern edge of 
the roughly-elliptical superbubble is twice that of the 
larger superbubble. If we ignore N186 D and concen- 
trate solely on the superbubble we see that the detected 
X-ray emission in the E — 0.3 - 2.8 keV energy band is 
relatively uniformly distributed, although with a slight 
tendency toward being brighter on the south and east- 
ern interior edges of the optical shell. 

The Hi column density along the line of sight to DEM 
L50 peaks in the projected center of the bubble, however 
(see Fig. |2fc), raising the possibility that the weak limb 
brightening is due to central absorption rather than true 
enhancement of the X-ray emission at the shell walls. 

That at least some of the Hi lies between us and the 
X-ray emitting plasma is clear from the way the X-ray- 
bright southeast limb of the superbubble fills in the gap 
in Hi column density at that location (compare Fig. [5^ 
and c). Based o n the observed velocities of the Hi gas, 
lOev et all 1|2002|) show that the Hi associated with DEM 
L50 is most likely the outer layer of the expanding shell 
of the bubble. 

DEM L50 is most prominent in the E = 0.6 - 1.6 
keV band. The SNR N186 D is prominent in both this 
energy band and the softer E = 0.3 - 0.6 keV energy 
band, despite being associated with a column density of 
Nh ~ 10 21 cm" 2 of Hi. This suggests that the SNR may 
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Fig. 3. — Adaptively-smoothed background-subtracted images of the diffuse X-ray emission associated with DEM L152. Panel (a) shows 
the E = 0.3 - 2.8 keV emission on a square root intensity scale, overlaid with contours increasing in factors of two beginning at a surface 
brightness of T,x = 4 X 10 — 7 c ounts s —1 arcsec -2 . These contour s are shown overlaid on a narrow-band Ha image fr om the Magellani c 
Emission Line Survey (MCELS: [Smith fc The MCELS Teamlll999f ) in (b). Panel (c) shows the Hi column density from lKim etU] (2003), 
linearly spaced between JVjj = 2 X 10 21 cm -2 (white) and 6 X 10 21 cm -2 (black), overlaid with the E = 0.3 - 2.8 keV emission contours 
from (a). Panels (d), (e) and (f) show adaptively smoothed X-ray images in the E = 0.3 - 0.6, E = 0.6 - 1.6 and E = 1.6 - 2.8 keV 
energy bands respectively, using a square root intensity scale. Contour levels increase in factors of two beginning at a surface brightness of 
T,x = 2 X 10~ 7 counts s _1 arcsec~ 2 . Dashed lines denote the edges of the Chandra ACIS CCD chips, and dashed red lines in (a) and (f) 
show the spectral extraction regions. 



be spectrally distinct from the superbubble. Neither the 
superbubble nor the SNR are apparent in the E — 1.6 - 
2.8 keV energy band. 

3.2.1. DEM L152 

In contrast to the relative simplicity of DEM L50, the 
spatial structure of DEM L152 (N44) in the X-ray, op- 
tical and Hi is complex (see Fig. 2]). The field of view 
of the Chandra observations (Fig. [3]) cover the following 
features id e ntifie d in the ROSAT PSPC observations of 
IChu et all (fl993T) : Shell 1, the main superbubble, and 
parts of Shell 3, a SNR, and the South Bar, in particu- 
lar the suspected break-out region where Shell 1 may be 
venting its contents through a break in its southern edge. 

The relationship between the new Chandra observa- 
tion and the features identified in prior X-ray observa- 
tions is shown in Fig. [5j where the Chandra ACIS field 
of view and contours of diffuse X-ray surface brightness 
are overlaid on a smoothed ROSAT HRl image of N44. 
The image was taken from the 108.5 ks HRI observation 
RH600913, obtained from the HEASARC data archive. 

Shell 1 is the most prominent X-ray source in the E = 
0.3 - 2.8 keV energy band Chandra image, and the soft 
diffuse X-ray emission is particularly bright around the 
inner edge of the optical shell. This limb-brightening is 
not an artifact of adaptive smoothing, as it is also clearly 
visible in images smoothed with a uniform FWHM= 30" 
Gaussian mask. 

The limb-to-center brightness contrast is a factor of 



~ 2 in the E = 0.3 - 2.8 keV and E = 0.6 - 1.6 keV 
energy band images and is slightly larger in the softer 
E — 0.3 - 0.6 keV energy band image (Fig. [3H). 

This central decrement in X-ray surface brightness can- 
not be due to absorption by intervening gas as both the 
X-ray-bright limbs and the fainter center of Shell 1 lie 
within a region of lower Hi column density (Fig. . 

Shell 3 falls partly within the ACIS S3 and S2 chips, 
although its eastern edge is not covered in the current 
observations. It has a comparable X-ray surface bright- 
ness to Shell 1 , even though the Hi column density along 
this line of sight (N-r ~5x 10 21 cm~ 2 ) is approximately 
twice the value toward the center of Shell 1. 

Bright, diffuse, soft X-ray emission also extends to 
both the west and south of S hell 1; the latter feature 
being known as the South Bar (|Chu et al.lll993f ). 

It is also noteworthy that both Shell 1 and Shell 3 
are clearly detected in the E — 1.6 - 2.8 keV energy 
band, in contrast to the complete lack of emission in this 
energy band from DEM L50 and N186 D. The South Bar 
lacks appreciable emission in this energy band, wh ich is 
consistent with the finding of iMagnier et al.l (|1996l) that 
the South Bar was spectrally softer than Shell 1. 



4. SPECTRAL ANALYSIS 



4.1. 



Spectral region selection 

As a first step to selecting regions for more detailed 
spectral analysis we constructed three-color images of the 
two bubbles where the resultant color is representative of 
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Fig. 4. — (a) A three-color optical/X-ray composite image of DEM L50 using narrow-band optical Ha+[Nn] emission (red), [Qu i] emission 
(green ) and diffuse soft X-ray emission in the E = 0.3-2.8 keV energy band (blue). The Ha, [Nil], and [Om] image comes from lOev et al.l 
(2002). The intensity in each band is shown on a square-root intensity scale, (b) A three-color soft X-ray composite image of DEM L50 
using the E = 0.3-0.6 keV (red), E = 0.6-1.6 keV (green) and E = 1.6-2.8 keV (blue) energy bands. The intensity in each band is shown 
on a linear intensity scale, (c) The statistical significance (in units of <r) of the dev iation of the local X-ray spectral hardness ratio, assuming 
uniform absorption, from the mean bubble hardness ratio, as described in § 14.11 overlaid with the contours of E = 0.3 — 2.8 keV diffuse 
X-ray surface brightness. Panels (d), (e) and (f) are the equ i valent of panels (a), (b) and (c), except that they show DEM L152. The 
Ha, [Nil] , and [Om] data is from Smith & The MCELS Team (1999). The field of view in all panels is the same 13' X 13' region shown in 
Figs. [2] and [3] 



TABLE 2 

Count Rates 



Target 


Region 


Area (arcmin 2 ) 


Rtot 


RdiR 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


DEM L50 


Bubble 


46.9443 


0.1830 ±0.0031 


0.1577 ±0.0030 


0.0252 ± 0.0009 


DEM L50 


SNR 


9.8197 


0.0776 ± 0.0018 


0.0749 ±0.0017 


0.0027 ± 0.0003 


DEM L152 


Bubble 


13.2994 


0.2633 ± 0.0041 


0.2626 ± 0.0041 


0.0006 ± 0.0004 


DEM L152 


SNR 


4.9459 


0.0422 ±0.0018 


0.0422 ± 0.0018 


0.0000 ±0.0001 


DEM L152 


South Bar 


11.8618 


0.1400 ±0.0031 


0.1362 ±0.0031 


0.0038 ± 0.0005 


DEM L152 


West 


17.4595 


0.1888 ±0.0037 


0.1855 ±0.0036 


0.0033 ± 0.0005 



the mean energy of the diffuse X-ray emission (Fig. 
& e) within the soft X-ray band. In these panels red 
represents the lowest energy photons (E = 0.3 - 0.6 keV), 
green intermediate energy photons (E = 0.6 - 1.6 keV) 
and blue relatively harder emission in the E = 1.6 - 2.8 
keV band. 

As we shall later show, these figures somewhat exag- 
gerate the spectral differences between the two objects 
and the degree of small-scale spectral variation within 
each object. Nevertheless they do serve some purpose by 



more clearly illustrating some of the possible spatial vari- 
ation in the distribution of diffuse emission at different 
energies and help inform our choice of regions over which 
to study the spectral properties of the diffuse emission. 

In Fig. [4p the emission from the superbubble DEM L50 
appears slightly spectrally harder than the brighter emis- 
sion from the SNR N186 D, but there is only a marginal 
suggestion of spectral hardness variation within DEM 
L50 itself. Based on this and the optical morphology of 
the bubble, we choose to separate the X-ray data into 
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Fig. 5. — A comparison between the the ROSAT HRI (grey, 
square root intensity scale) and Chandra ACIS observations (con- 
tours, as in Fig. [3^,) of N44. Dashed lines denote the e dges of the 
Chand ra ACIS CCD chips. The features identified in IChu et aXI 
(1993) are outlined and labeled. 

two regions for spectral fitting: An elliptical region en- 
compassing SNR N186 D and another elliptical region 
covering the superbubble but excluding the SNR. These 
are the regions shown in Fig. [5^. In order to examine 
spectral variations on smaller scales, we further subdi- 
vided DEM L50 into a Center and Limb region when 
performing spectral fits (see § 14. 2[) . Note that only the 
parts of those ellipses within the S3 chip are actually 
included within the spectral extraction regions. 

The soft X-ray emission in the vicinity of DEM L152 
shows more evident spectral variation than does DEM 
L50 (Figure Sb). Based on this and the optical mor- 
phology of the N44 nebula we chose s pectral extraction 
regi ons similar to those us ed earlier bv lChu et al.l (|1993| ) 
and iMagnier et all (|1996Q : Shell 1 (the main bubble), 
Shell 3 (the SNR), and the South Bar, with most of the 
remaining emission associated with a new region we sim- 
ply name West (see Fig. [3k). 

Diffuse X-ray count rates within the spectral extraction 
regions of both sets of observations are given in Table [2j 
Column 3 gives the area of each region. Columns 4-6 
show the background-subtracted ACIS S3 count rates in 
the E = 0.3 — 2.8 keV energy band within each specified 
region, where i? to t gives the count rate for the full X- 
ray emission, i?diff gives the estimated diffuse emission 
count rate, and i?ps gives the total point source count 
rate. Emission from point-like X-ray sources accounts 
for ~ 14% of the total soft X-ray emission within the 
boundaries of the DEM L50 bubble, and less than 1% 
of the emission within the boundary of the DEM L152 
Shell 1. 

Given Chandra's high spatial resolution and sensitiv- 
ity we should be able to determine whether the spectral 
properties of the diffuse emission vary on spatial scales 
smaller than the entire shell or bubbles. In Figures 0t> 
& e we show energy-color-coded images of the diffuse 
emission around DEM L50 and DEM L152. Such im- 
ages can be prone to smoothing-related artifacts, so we 
also mapped the hardness ratios in I' x 1' pixels. These 



large pixels were required in order to minimize the sta- 
tistical uncertainty in the hardness ratio in each pixel. 
The spectral hardness ratio, which will be discussed in 
more detail in §[00 is defined as Q = (H - S)/(H+ S), 
where H and S are the count rates in the specified hard 
(E = 0.6 - 1.6 keV) and soft {E = 0.3 - 0.6 keV) energy 
bands respectively. Thus higher values of Q correspond 
to a harder spectrum. Images of the hardness ratio itself 
tend to be difficult to interpret, as they can be dominated 
by noise in the lower signal-to-noise regions. For this rea- 
son we constructed maps of the statistical significance of 
the deviation of the local hardness ratio from the average 
hardness ratio Q over each bubble. For a given pixel with 
hardness ratio Q P i X and uncertainty <7Q jP i x , Figures^]: & 

f show (Qpi x — Q)/cQ,pix- In other words, the value in 
each pixel corresponds to the significance of the devia- 
tion in Gaussian a. While Figures |4b & f hint at some 
of the intrinsic large-scale variations across the super- 
bubbles, variations in column density and the resultant 
absorption of soft X-ray photons have a substantial effect 
on the observed hardness ratios. Figures Hh & f should 
only be used to identify regions that appear spectrally 
distinct from one another. 

Although some spatial variation in the spectral prop- 
erties of the diffuse X-ray emission is present in both 
objects, it is not particularly strong, and emission in the 
E = 0.6 - 1.6 keV energy band dominates both sets of 
observations. Hardness ratio variations are present on 
large scales, for example between the bubbles and the 
nearby supernova remnant, or between Shell 1 and the 
Southern Bar region of DEM L152 . On smaller angu- 
lar scales (~ 1') we find that any spatial variations in 
spectral hardness are of marginal statistical significance 
in our current X-ray data. 

4.2. Spectral fitting 

For the spectral fits, we removed background flares 
from both the S3 and S2 chips with the LC_CLEAN routine 
in CIAO. We then used the re- normalized ACIS stowed 
background files to subtract the approximately constant 
particle background from both chips. To re-normalize 
the background files, we compared the relative count 
rates of the stowed background and our observations be- 
tween 8 keV and 9.4 keV and re-normalized the particle 
background to match the count rate of our data. The 
8 — 9.4 keV energy range was chosen to exclude real X- 
ray photon events and a line feature above 9.4 keV. 

We used the SPECEXTRACT script in CIAO to extract 
spectra from the specified regions and subtract the par- 
ticle background. To account for the X-ray background, 
we fit the extracted s p ectra from the S2 chip, follow- 
ing iKuntz fc Snowdenl ()2000l ) and iHenlev et all (|2007t ). 
and used the resulting fits as a component of our subse- 
quent models to the S3 data. A thermal plasma in col- 
lisional ionization equilibrium fits the local X-ray back- 
ground, while two additional thermal plasmas in equi- 
librium are used for the Galactic background, and the 
cumulative emission from extragalactic sources is fit by 
a power law. Our best fits to the S2 background are 
shown in Table |3J The power law index was fixed to the 
value found bv lChen et al.l (|1997| ). and th e column densi- 
ties w ere obtained from HEASARC from Kalbcrla et al. 
( 2005). Temperatures were left as a free parameter. The 
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Local Hot Bubble temperature is given in Column 3, 
the Galactic Hi column density is shown in Column 4, 
the temperatures of the two thermal components of the 
Galactic background are given in Columns 5 and 6, and 
the extragalactic power law index is shown in Column 7. 
Normalizations were re-scaled to account for the relative 
difference in area between the S3 extracted regions and 
the S2 chip. For the sake of comparison, we include two 
fits to the S2 background for DEM L50; the "Low S2" fit 
was a local minimum in x 2 before the global best fit was 
found. The effect of varying the background fit can be 
seen by comparing the DEM L50 models with the "High 
S2" and "Low S2" background fits. 

All extracted regions are shown in Figure [6] In ad- 
dition to the regions in Figs. [2^ and [3^, we extracted 
spectra for smaller subregions. We divided DEM L50 
into Limb and Center regions a separate spectrum for the 
bright southern portion of the Limb. For DEM L152, we 
extracted additional spectra for the Bubble (that is, Shell 
1) Limb, West Limb, and West Blowout regions as well 
as for the two knots shown in Fig.[6l While some clumps 
seen in the smoothed image of DEM L152 are artifacts 
of the smoothing algorithm, these two knots are statisti- 
cally significant (see Fig. 0f ) and appear when different 
smoothing algorithms are used. 

We fit each superbubble region and subregion with 
one- and two-temperature plasmas in collisional ioniza- 
tion equilibrium, using the XSPE C Astrophysical P lasma 
Emission Code (APEC) models (Sm ith et al.ll2001l ) with 
variable abundances (Tables 0] and [6]). We examined an 
additional, non-equilibrium ionization model for DEM 
L50 (Table 0. In all models, He, C, N, Al, and Ni abun- 
dances were fixed to the solar values. We varied the a el- 
ements O, Ne, Mg, Si, S, Ar, and Ca in proportion to one 
another, and Fe was left as a free parameter. Although 
absolute abundances are not well constrained, the a/Fe 
ratio is more reliable. To ensure that the fits were not 
local minima, we looked at possible values of Ah, tem- 
perature, and abundances spanning one to two orders of 
magnitude with the XSPEC command STEPPAR. As a 
further check on the validity of our fits and to establish 
error bars for each variable, we ran the error command 
on the free parameters. 

To determine fluxes, we ran the XSPEC command 
flux on each model and determined the absorption- 
corrected fluxes by setting Ah to zero. We adopted 
the distance of 50 kpc to the LMC when converting 
absorption-corrected fluxes to luminosities. Observed 
fluxes and the calculated luminosities are shown in Ta- 
bles [7]fTUl which we discuss below. 

We also calculated hardness ratios and median ener- 
gies for the diffuse emission from each region, shown in 
Table Q3J the hardness ratio Q = (H - S)/{H + S), 
where H and S are the luminosities in the hard (E = 0.6 
- 1.6 keV) energy band and soft (E = 0.3 - 0.6 keV) 
energy band. Absorption by intervening Hi results in 
an apparently harder spectrum, and Figs. [2]: & [3J: show 
that there is significant variation in the Hi column den- 
sity across these objects. The column density was a 
free parameter in our fits, and the resulting luminosi- 
ties and hardness ratios account for Ah variations be- 
tween regions. Q obtained using the count rate (i.e. not 
absorption-corrected) is shown for comparison in Column 
4 of Table [TT1 The influence of the column density can be 



seen by comparing the luminosity-determined and count 
rate-determined hardness ratios for the SNR 0523-679 
(jChu et al.lll99l in DEM L152. While the SNR appears 
to have the hardest spectrum from the count rates, its 
higher absorption column density masks a much softer 
s pectrum. 

iHong et al. (2004) demonstrate that the median energy 
is a more reliable indicator of X-ray spectral variations 
than the hardness ratio is. As a result, we include the 
median energy of each region in Table 111! The median 
energies generally agree with the hardness ratios. An 
exception is the DEM L152 SNR, which has a higher 
median energy due to the Hi absorption discussed above. 

4.3. DEM L50 

Modeled fits to DEM L50 are shown in Tables [4] and EJ 
The fits in Table |4] model the emission as one- and two- 
temperature thermal plasmas. The best fit Ah in Col- 
umn 3 is con sistent with the o bservations of the Hi in 
DEM L50 bv lOev et all (|2002D . The first temperature 
component, fcTi, appears in Column 4, with its normal- 
ization in Column 5, and the second component, kTi, 
appears in Column 6, with its normalization in Column 
7. Column 8 shows the "Effective kT", defined as the 
average kT of the region weighted by the normalization. 
The "Effective kT" thus shows the dominant tempera- 
ture component in the region. Columns 9 and 10 give 
the abundances relative to solar, with the a/Fe ratio in 
Column 11. Columns 12 and 13 show the reduced % 2 and 
the number of degrees of freedom of the fit. The statisti- 
cally best fits, determined by running ftest in XSPEC, 
are shown in bold. Figure [7] shows a sample fit to the 
Center region of DEM L50. 

From the best fits for the Limb and Center regions 
in Table |H the Center's temperature distribution ap- 
pears more homogeneous than the Limb's, and the two- 
temperature best-fit to the Limb shows a higher tem- 
perature component that may correspond to the Bright 
Limb subregion. This correlation would be consistent 
with shock-heating of the southern part of the Limb by 
an expanding SNR, while the rest of the Limb remains 
at a cooler temperature. 

4.3.1. Low a/Fe vs. non- equilibrium ionization 

In each region of DEM L50, spectral fits consistently 
gave low a/Fe ratios (see Tabled]), with typical values 
of about 0.3 compared to an average LMC value of 0.6 
(|Smithlll999t ). These low ratios are especially surprising, 
given that we might have expected a high a/Fe ratio due 
to an enhancement in a elements from core- collapse SNe . 
The bubble has likely experienced ~ 2 SNe (|Oevlll996bl) . 
so an a enhancement would be natural. Raising the a/Fe 
ratio by forcing the abundances to solar values resulted 
in the models labeled 'solar' in Table |4] In every case, 
ftest showed that the free abundance fits were statisti- 
cally better than the solar values. The ftest command 
calculates a p-value, which here denotes the probabil- 
ity that solar-abundance material could randomly pro- 
duce a spectrum fit by the more complicated, free abun- 
dance fit. The resulting p-values ranged from 0.049 to 
as low as 4 x 10 -6 , showing that the solar fits have a 
very low probability of being correct. Although the solar 
abundances are statistically excluded, there may be sys- 
tematic problems in forward-fitting low-resolution X-ray 
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TABLE 3 

Parameters for Background Fits 



Bubble 


Background Fit 


fcT LHB (keV) 


N H a (10 22 cm" 2 ) 


feTi (keV) 


kT 2 (keV) 


Photon Index b 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


DEM L50 


High S2 Fit 


1.08 


0.132 


0.05 


0.22 


1.46 


DEM L50 


Low S2 Fit 


0.08 


0.132 


0.10 


0.24 


1.46 


DEM L152 


L152 S2 Fit 


0.08 


0.228 


0.06 


0.25 


1.46 



a From lKaTberiTc t al. (200^), provided by HEASARC. 
b From lChen et al.l U997I) . 

spectra, which could give us misleading estimates of the 
abundances and their uncertainties. Using the Mckal or 
Raymond-Smith models instead of APEC did not affect 
the a/Fe ratios. 

Another possibility is that collisional ionization equi- 
librium is not a valid assumption for DEM L50. To test 
this scenario, we performed the non-equilibrium ioniza- 
tion fits shown in Table [5] Temperature is shown in 
Column 4, the normalization is in Column 5, abundances 
relative to solar are in Columns 8 and 9, and the reduced 
X 2 and number of degrees of freedom are in Columns 11 
and 12, respectively. The a/Fe ratios found with these 
new fits are all consistent with solar values. A sample 
non-equilibrium ionization fit is shown in Figure [8] for 
the DEM L50 Center region. 

To see whether or not non-ionization equilibrium is a 
plausible scenario, we consider the ionization timescale. 
The ionization timescale in seconds is found by dividing 
r, the density-weighted ionization timescale in s cm~ 3 
(Column 6), by the electron density, n e . To find the 
electron density, we used 

10~ 14 f 

where r\ is the normalization of the fit, in units of cm' 5 
and Da is the distance to the superbubble, assumed 
to be 50 kpc. We also assumed a spherical geometry, 
uniform density distribution, and ^10% He abundance 



(n e ~ 1.2 hh; see, e.g.. lPeimbert et aT1l2007|) . The result- 
ing ionization timescale estimates are shown in Column 
7 of Tableland are all on the order of 10 5 years, about 
10% of the superbubble's age. While our observation of 
the bubble in this stage would be a chance occurrence 
to some extent, this timescale is not unrealistically brief 
and a recent supernova is also consistent with many of 
the observed properties of the bubble as we will discuss 
later. Furthermore, the hypothesized brightening from 
SNR impacts should enhance the detectability of bubbles 
with recent SNe, increasing the likelihood of observing a 
bubble that is not in equilibrium. 

If the low a/Fe abundances are indeed real, they may 
reflect an Fe enrichment in the ISM near DEM L50. 
DEM L50 i s located in a l arge, 1.4 kpc diameter void 
in the ISM (|Oev et al.ll2002f ) ; a single Type la supernova 
could not noticeably enrich this entire region, nor could 
it enrich the esti mated 1.7 x 10 5 Mq of material swept 
up by DEM L50 (|Oev et al.ll2"002l) . On the other hand, if 
the density is low enough, the interior of DEM L50 could 
be enriched by a single Type la event, and SNR N186 D 
near the bubble would be a likely candidate. The a/Fe 
ratio of the SNR, however, does not seem noticeably dif- 
ferent from the rest of the bubble (see § 14. 5p and may in 
fact be less Fe-enriched than the Center region. There is 
no clear spectral contrast between the SNR and regions 
dominated by the shocked ISM like tha t seen in the Type 
la SNR DEM L71 (|Hughes et al.ll2003| ). nor are there any 
noticeable Fe, S, or Si features such as those seen in Type 



TABLE 4 
Fit Parameters for DEM L50 



Model 


Bkgd Fit a 


Af H (10 22 cm- 2 ) 


fcTi (kcV) 


normi 


kT 2 (kcV) 


norm2 


eff kT (keV) 


a/O© 


Fe/Fe 


a/Fe 
(a/Fen) 
(11) 


red. x 2 


D.O.F 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(12) 


(13) 


Center IT 


High S2 


n osj + 0.03 
UcS8 -0.04 


u±£> -0.01 


3.5E-03 






0.15 


0.88t°;! 7 




0.l5lS;g» 


1.743 


211 


Center IT solar 


High S2 


0.31 ± 0.05 


0.19 ± 0.01 


1.3E-03 






0.19 


1.00 


1.00 


1.00 


1.943 


213 


Center 2T 


High S2 


40+ ' 08 
u.+u_o Qg 


u - li -o.n 


4.7E-03 


0.17 ± 0.03 


1.6E-03 


0.13 


1 10+ 1 ' 18 


6.03t 3 ° 3 f 


n 18+°- 14 


1.757 


209 


Center 2T solar 


High S2 


n o 9 +0.06 


ll+ - 04 

u - 11 _0.02 


3.1E-03 


n 9 o+0.06 
u - zo -0.02 


5.9E-04 


0.13 


1.00 


1.00 


1.00 


1.835 


211 


Limb IT 


High S2 


0.29 ± 0.05 


n 1s +0.02 
0.1S_ .04 


2.7E-03 






0.18 


54+ - 29 


1 a 3 +8.04 


0-28l«; ° 


1.744 


302 


Limb 2T 


High S2 


n ,, + 0.03 


n ns + 0.20 
UU8 -0.08 


1.4E-02 


n 22+° 07 


7.3E-04 


0.09 


-. „_+0.62 
1 - ao -0.36 


2.7811;^ 


0.45t8# 


1.718 


300 


Limb 2T solar 


High S2 


0.21 ± 0.04 


11+ - 04 


1.4E-03 


n 9 7 +0.03 
u - z '-0.04 


2.1E-04 


0.13 


1.00 


1.00 


1.00 


1.749 


302 


Limb IT 


Low S2 


0.36 ± 0.06 


U - 18 -0.03 


2.6E-03 






0.18 


n QC -+2.70 
0-95_ 39 


O ,0 + 13.92 

o.oo i g$ 


n 90+0.12 


1.793 


302 


Limb 2T 


Low S2 


o.34t°;?I 


14+ - 26 
u - ±4 -0.14 


1.4E-03 


21 + ' 16 
u - zl -0.11 


7.1E-04 


0.16 


-1 o 7 +2.50 
± - J '-0.71 


3.58l»;» 


0-38l°; 2 « 


1.799 


300 


Bright Limb IT 


High S2 




n 1Q+0 O 2 
ul3 -0.03 


1.9E-03 






0.19 


o.56tg:SS 


1.07111° 


n -o+0. 27 


1.516 


137 


Bright Limb 2T 


High S2 


0.37 ± 0.07 


n 11+0.16 
U - li -().ll 


1.8E-03 


22+ ' 11 
u - zz -0.05 


8.3E-04 


0.14 


o-8o±S:« 


, 91 +2.72 

1 - zl -0.72 


cc+0.44 

0.66I 037 


1.532 


135 


SNR 2T 


High S2 


0.09 ±0.02 


0.19 ±0.01 


4.7E-04 


n -9 + 0.09 
' "^-0.05 


2.6E-05 


0.22 


0.37l° ; 2 ? 


■1 C(3 + 2.21 


n 22+ 13 


2.114 


145 


SNR 2T solar 


High S2 


20+ 06 
u - zu -0.03 


0.08±g;S 


6.2E-03 


0.29 ± 0.02 


2.8E-04 


0.09 


1.00 


1.00 


1.00 


2.173 


147 


SNR 2T 


Low S2 


06+ 05 
u - UD -0.04 


0.19±°;°? 


4.3E-04 


0.73 ± 0.09 


3.2E-05 


0.23 


27+ - 12 

u - z '-0.06 


1 26+ 1 ' 50 


22 +o n 
u - zz -0.06 


2.150 


145 


SNR 2T solar 


Low S2 


25+ - 04 


0-08±g;g| 


9.9E-03 


27+ ' 02 
u -^'-o.oi 


3.9E-04 


0.09 


1.00 


1.00 


1.00 


2.255 


148 



CO 

p 

a 
*-i 

P 

cr 
aT 

P 

c 

TO 

Cd 
p 





Note. — Errors show the 90% confidence interval when one parameter is varied. 
a The Low S2 background gave poor fits to the Center region which are not included. 



TABLE 5 
DEM L50 Non-Equilibrium Fits 



Model 
(1) 



Bkgd Fit a JV H (10 22 cm- 2 ) fcTi (kcV) normi 
(2) (3) (4) (5) 



t (s cm 3 ) 
(6) 



Ion. Time, (yr) 
(7) 



a/0 Q 
(8) 



Fe/Fe 
(9) 



a/Fe 
(a/Fe Q ) 
(10) 



red. x 2 D.O.F. 



(11) 



(12) 



Center NEI High S2 

Limb NEI High S2 

Limb NEI Low S2 

Bright Limb NEI High S2 

SNR NEI High S2 



0.3V 



1+0.05 
-0.08 
n ok+0.04 
°- 26 -0.03 



0.34" 1 



^+0.09 
-0.08 

34+ - 12 



0.11 



+0.05 
-0.03 



0.28 
0.34 
0.28 



0.07 
0.06 
+0.11 
-0.05 
+0.11 
0.07 



S2+ ' 14 

U-^-O.lO 



0.58 



+0.23 



5.r 

3.2 
5.5 
3.8 



J E-04 5.57E10 ± 4.68E10 22,000 



E-04 



-2.3 

! 4 7 E-04 



+ 18.7 

4 



3.31E10 ± 2.84E10 46,000 
6.26E10 ± 4.69E10 65,000 



+15 ' 8 E-04 8.33E10 ± 2.88E10 53,000 b 

i+0.04 T 



1.2I-.-E-04 



2.60E10 ± 6.39E9 11,000 



76+ 36 

U -' D _0.24 



0.83 
1.28 
0.70 
0.29 



+0.46 
-0.24 
+ 1.48 



+0.69 
-0.36 
+0.10 



74+0-36 
u -' 4 -0.26 

1 in - 62 

i - lu -0.34 



1.511 
0.58: 
0.45 



0.60 

0.55 



+0.14 



1.043 
0.76 
0.85 
1.21 



0.643;° o: « 



+ 0.16 
-0.12 
+0.19 



+0.38 



1.838 
1.718 
1.718 
1.543 
2.265 



210 
301 
301 
136 
146 



Note. — Errors show the 90% confidence interval when one parameter is varied, r errors are estimates only. 
a The Low S2 background gave poor fits to the Center region which are not included. 
b assuming n e ~ 0.05 cm -3 . 



TABLE 6 
Fit Parameters for DEM L152 



Region and Model 


at /io22 — 2\ 

Nn(10 cm z ) 


kTi (keV) 


normi 


kT2 (keV) 


norni2 


eff kT (keV) 


«/o 


Fe/Fe Q 


a/Fe 
(a/Fe ) 


j 2 

red. x 


D.O.F 








(A\ 
\V 




W 












Bubble IT 


0.27 ± 0.02 


0.37 ±0.02 


2.5E-03 






0.37 


o.68tS;IS 


„ on +0.16 
UJ1 -0.07 


o o7 -+0.58 
J Z ' -0.43 


1.360 


151 


Bubble 2T 


29+ ' 04 


17+ - 14 
u - 1 '-0.07 


5.5E-04 


n o 7 +0.05 
— 0.02 


2.2E-03 


0.33 


0.76tg;™ 


26+ ' 12 


2 95+ ' 65 
Z - MO -0.47 


1.352 


149 


Bubble Limb IT 


0.28 ±0.03 


0.33 ±0.02 


3.3E-03 






0.33 


0.48i°;?? 


0.18 ±0.05 


^• OO -0.43 


1.794 


138 


Bubble Limb 2T N H fixed 


0.28 


14+"- 10 


4.7E-04 


S5+ ' 02 


1.7E-03 


0.30 




n o 9 +0.39 
U - Jz -0.15 


o 7^+0.55 
z -'°-0.38 


1.782 


137 


West IT 


n 97 +0.03 
u - z '-0.02 


0.30 ± 0.02 


3.6E-03 






0.30 


0-36t°;? 4 


o in+ 007 


3 42+ 1 ' 34 


1.088 


154 


West 2T 


n 9s +0.04 


UZO -0.03 


2.5E-03 


0.79l° - 


2.9E-04 


0.31 


0.60tg;« 


29+° 29 


1 74+0 97 
A - ' -0.58 


1.027 


152 


West Limb 1T C 


n 24+ - 06 

u - z4 -0.03 


S2+ 04 


7.3E-04 






0.32 


u -+ y -0.29 


n ni +2.28 
u - zi -0.09 


o ,1+1.21 
z - 01 -1.89 


1.387 


78 


West Limb 2T N H fixcd c 


0.28 


21+ 03 


6.6E-04 


0.83in 9 4 


5.7E-05 


0.26 




i r>o+2.33 
1 - zz -0.75 


"S6+ - 35 
u - aD -0.22 


1.355 


77 


West Blowout IT N H fixed 


0.28 


0.25 ±0.02 


2.0E-04 






0.25 


2.34 b 


U,iZ -0.32 


7.39-4.22 


1.143 


77 


West Blowout 2T N H fixed 


0.28 


14+ 4 - 30 


2.5E-05 




3.7E-05 


0.22 


10.10 b 


1 6 4+374.11 


6.08-3.23 


1.156 


75 


South Bar 1T C 


0.23 ± 0.02 


n otr+0.03 
u - od -0.02 


1.0E-03 






0.35 


75+ 2 - 71 
u - '°-0.35 


21 + ' 70 
u - zl 0.08 


3 56+ 1 ' 35 


1.054 


128 


South Bar 2T C 


0.21± 0.02 


29+ 04 


4.4E-05 




1.5E-05 


0.43 


lO.OOt 21 ^ 2 


6.84 b 




1.034 


126 


SNR IT 


n 34+0 06 


n „ s +0.07 
Ud8 -0.06 


9.3E-04 






0.38 


0.36t™? 


n 04+° 05 
uu *-0.04 


10.24+^f 


1.176 


67 


SNR 2T 


U -J°_0.04 


n o.r+0.08 
U.OO_ Q 23 


7.4E-04 


0.86 b 


6.3E-05 


0.39 


0.44±°;£ 


o m +om 

U - UO -0.05 


o 1H +51.18 


1.198 


65 


Knot 1 IT solar 


0.22 


0.30 ± 0.03 


8.2E-05 






0.30 


1.00 


1.00 


1.00 


1.471 


17 


Knot 1 IT free abun 


0.22 


0.30 ± 0.03 


7.3E-05 






0.30 


1.40-1. 07 b 


5S+ 9 ' 87 


2 fi3+ 4 ' 67 


1.253 


15 


Knot 2 IT solar 


0.22 


S0+ 07 


3.5E-05 






0.30 


1.00 


1.00 


1.00 


1.288 


8 



CO 
P 
"O 
o 
>-j 
cr 
pi 

cr 
cr 
ST 

■ 

P 

CD 
TO 

Cd 
p 
a 

aq 



Note. — Errors show the 90% confidence interval when one parameter is varied. 
b Errors arc unconstrained. 

c It is not clear which fit is statistically better. The p-value for a 2T fit is ~ 0.1. 
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TABLE 7 
Observed Fluxes for DEM L50 



Model 
(1) 


0.3-2.0 keV 
(2) 


0.5-2.0 keV 

(3) 


0.3-8.0 keV 
(4) 


0.5-8.0 keV 

(5) 


2.0-8.0 keV 

(6) 


Center IT 


0.24 


0.23 


0.28 


0.26 


0.032 


Center IT solar 


0.24 


0.22 


0.27 


0.25 


0.032 


Center 2T 


0.24 


0.23 


0.28 


0.26 


0.031 


Center 2T solar 


0.25 


0.23 


0.28 


0.26 


0.032 


Center NEI 


0.25 


0.22 


0.28 


0.26 


0.032 


Limb IT 


0.39 


0.36 


0.45 


0.42 


0.057 


T imh 9T 


0.41 


0.36 


0.47 


0.42 


0.057 


Limb 2T solar 


0.43 


0.36 


0.48 


0.42 


0.058 


Limb IT (LowS2) 


0.48 


0.36 


0.54 


0.43 


0.063 


Limb 2T (LowS2) 


0.55 


0.36 


0.55 


0.43 


0.064 


Limb NEI 


0.41 


0.36 


0.46 


0.42 


0.057 


Limb NEI (LowS2) 


0.48 


0.36 


0.55 


0.43 


0.064 


Bright Limb IT 


0.16 


0.15 


0.17 


0.17 


0.015 


Bright Limb 2T 


0.16 


0.15 


0.18 


0.17 


0.015 


Bright Limb NEI 


0.16 


0.15 


0.18 


0.17 


0.015 


SNR 2T 


0.37 


0.26 


0.39 


0.28 


0.021 


SNR 2T solar 


0.34 


0.26 


0.36 


0.28 


0.020 


SNR 2T (LowS2) 


0.36 


0.26 


0.39 


0.29 


0.023 


SNR 2T solar (LowS2) 


0.35 


0.26 


0.37 


0.28 


0.022 


SNR NEI 


0.35 


0.26 


0.37 


0.28 


0.021 


Note. — All fluxes arc in units of 


10~ 12 erg cm" 


- 2 s _1 . LowS2 


indicates that 


the Low S2 



background fit was used. 



TABLE 8 

Luminosities and Surface Brightnesses for DEM L50 



Model 


0.3-2.0 keV 


0.5-2.0 keV 


0.3-8.0 keV 


0.5-8.0 keV 


2.0-8.0 keV 


Surface Brightness 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Center IT 


15.69 


8.43 


15.79 


8.53 


0.098 


2.70 


Center IT solar 


9.41 


4.93 


9.51 


5.03 


0.099 


1.63 


Center 2T 


21.02 


10.51 


21.12 


10.61 


0.098 


3.61 


Center 2T solar 


12.41 


5.47 


12.51 


5.57 


0.098 


2.14 


Center NEI 


10.25 


4.95 


10.35 


5.05 


0.099 


1.77 


Limb IT 


14.42 


6.67 


14.60 


6.85 


0.18 


1.40 


Limb 2T 


23.81 


8.25 


23.98 


8.42 


0.18 


2.31 


Limb 2T solar 


10.95 


3.99 


11.13 


4.16 


0.18 


1.07 


Limb IT (LowS2) 


13.72 


9.60 


13.92 


9.80 


0.20 


1.32 


Limb 2T (LowS2) 


12.54 


8.73 


12.74 


8.93 


0.20 


1.22 


Limb NEI 


12.54 


5.35 


12.72 


5.53 


0.18 


1.22 


Limb NEI (LowS2) 


12.15 


8.48 


12.35 


8.68 


0.20 


1.19 


Bright Limb IT 


7.72 


4.50 


7.77 


4.54 


0.047 


2.81 


Bright Limb 2T 


8.50 


4.33 


8.54 


4.38 


0.047 


3.08 


Bright Limb NEI 


5.83 


3.43 


5.87 


3.48 


0.048 


2.12 


SNR 2T 


3.45 


1.35 


3.51 


1.41 


0.063 


1.02 


SNR 2T solar 


8.69 


2.73 


8.75 


2.79 


0.060 


2.54 


SNR 2T (LowS2) 


1.70 


1.11 


1.77 


1.18 


0.070 


0.51 


SNR 2T solar (LowS2) 


10.83 


3.81 


10.90 


3.88 


0.067 


3.16 


SNR NEI 


3.67 


1.47 


3.73 


1.54 


0.064 


1.08 


Totals of IT Models 


30.11 


15.11 


30.39 


15.38 


0.27 


1.87 


Totals of 2T Models 


44.83 


18.76 


45.10 


19.03 


0.28 


2.78 


Totals of NEI models 


22.79 


10.30 


23.07 


10.58 


0.28 


1.42 



Note. — Luminosities are in units of 10 35 crgss -1 . Surface brightness is shown for the 0.3 — 8.0 keV band in units 
of 10 31 erg s _1 arcscc -2 . LowS2 indicates that the Low S2 background fit was used. Totals are sums of Center and Limb 
regions, with High S2 background fits, free Nu, and free abundances. 
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TABLE 9 

Observed Fluxes for DEM L152 



Region and Model 


0.3-2.0 kcV 


0.5-2.0 kcV 


0.3-8.0 kcV 


0.5-8.0 keV 


2.0-8.0 kcV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


T 1 1 1 1 1 ' 1 ' 

rJuDble ±1 


0.93 


0.87 


0.95 


0.89 


0.016 


rSubble 1 1 


0.93 


0.88 


0.94 


0.89 


0.017 


1.1.1.. T 1 T" 

rSubble L*imb 11 


0.77 


0.72 


0.78 


0.73 


n /www.' 

u.uuyo 


rSubble Limb zl ivjr nxea 


0.77 


0.73 


0.78 


0.74 


0.011 


West 1 1 


0.73 


0.67 


0.73 


0.67 


0.0061 


west A 1 


U. I o 


U.OO 


U.io 


U. (U 




West Limb IT 


0.26 


0.23 


0.26 


0.23 


0.0020 


West Limb 2T JV H fixed 


0.25 


0.24 


0.26 


0.24 


0.0061 


West Blowout IT jV H Axed 


0.19 


0.18 


0.20 


0.18 


0.00070 


West Blowout 2T AT H fixed 


0.20 


0.18 


0.20 


0.18 


0.00077 


South Bar IT 


0.53 


0.48 


0.54 


0.48 


0.0053 


South Bar 2T 


0.54 


0.48 


0.56 


0.49 


0.016 


SNR IT 


0.15 


0.14 


0.16 


0.15 


0.0042 


SNR 2T 


0.15 


0.15 


0.16 


0.15 


0.0076 


Knot 1 IT solar 


0.049 


0.046 


0.049 


0.046 


0.00019 


Knot 1 IT free abun 


0.051 


0.048 


0.052 


0.048 


0.00036 


Knot 2 IT 


0.024 


0.022 


0.024 


0.022 


0.00012 



Note. — Fluxes are in units of 10 12 erg cm 2 s 1 . 



TABLE 10 

Luminosities and Surface Brightnesses for DEM L152 



Region and Model 
(1) 


0.3-2.0 keV 
(2) 


0.5-2.0 keV 

(3) 


0.3-8.0 keV 
(4) 


0.5-8.0 kcV 
(5) 


2.0-8.0 kcV 
(6) 


Surface Brightness 
(7) 


Bubble IT 


20.20 


9.10 


20.26 


9.16 


0.054 


4.23 


Bubble 2T 


22.22 


10.58 


22.27 


10.64 


0.056 


4.65 


Bubble Limb IT 


17.88 


8.61 


17.91 


8.64 


0.032 


4.31 


Bubble Limb 2T iV H fixed 


17.66 


8.54 


17.69 


8.58 


0.035 


4.26 


West IT 


21.59 


8.50 


21.61 


8.52 


0.020 


3.44 


West 2T 


22.19 


9.12 


22.27 


9.20 


0.075 


3.55 


West Limb IT 


6.23 


2.42 


6.24 


2.43 


0.0066 


3.26 


West Limb 2T jV H fixed 


7.31 


3.24 


7.33 


3.26 


0.020 


3.83 


West Blowout IT JV H fixed 


7.22 


2.70 


7.22 


2.71 


0.0023 


3.07 


West Blowout 2T JV H fixed 


7.26 


2.79 


7.27 


2.79 


0.0025 


3.09 


South Bar IT 


12.81 


4.72 


12.82 


4.73 


0.017 


3.04 


South Bar 2T 


11.99 


4.29 


12.04 


4.34 


0.050 


2.85 


SNR IT 


5.72 


2.12 


5.73 


2.14 


0.014 


3.22 


SNR 2T 


5.72 


2.16 


5.75 


2.18 


0.025 


3.23 


Knot 1 IT solar 


0.81 


0.45 


0.81 


0.45 


0.0006 


3.70 


Knot 1 IT free abun 


0.82 


0.46 


0.82 


0.46 


0.0012 


3.74 


Knot 2 IT 


0.43 


0.21 


0.43 


0.21 


0.0004 


4.02 


Total of IT Models 


54.60 


22.32 


54.69 


22.41 


0.091 


3.58 


Total of 2T Models 


56.40 


24.03 


56.58 


24.21 


0.18 


3.70 



Note. — Luminosities are in units of 10 35 ergss 1 . Surface brightness is shown for the 0.3 — 8.0 keV band in units of 10 31 erg 
s _1 arcsec — 2 . Totals are sums of Bubble, West, and South Bar regions. 



Ia SNR N103B (IHuehes et al.l[l995lh Furthermore, SNR 
N186D's proximity to young stars and star-forming re- 
gions make a core-coll apse SN origin more likely (e.g. 
IChu fc Kennlxut^[T98l . 

Dust destruction from a supernova shock wave can- 
not explain the low a/Fe ratio. Although dust destruc- 
tion would increase the Fe abundance, the destruction 
of silicates would concomitantly raise the O ab undance 
(Compiegne, M., private communication; see iWhittetj 
2003 for O and Fe abundances in dust). We note that 
DEM L50 does have an uncommon infrared morphology 
(Slater et al. 2010, in prep.); its 8 fim emission appears 
to fill the bubble's Ha shell, whereas many other LMC 
superbubbles only show polycyclic aromatic hydrocarbon 
(PAH) emission external to the bubble. This morphol- 
ogy is likely explained by PAH emission on the interface 



between the ionized and the foreground neutral layer, 
which is seen in projection toward the center of the bub- 
ble (see Fig. [2] and § 13 . 2[) . An unusual dust environment 
does not appear to be the cause of the low observed a/Fe 
ratio. 

Similarly odd Fe abundances have been observed in 
other objects. DEM L316 consists of a pair of interacting 
SNRs; on e SNR has a high Fe abundance while the othe r 
does not (|Williams fc Chul [20051: iNishiuchi et al.1 12001ft . 
As mentioned previously, t he LMC a/Fe ratio is low in 
general, -0.6 (Smith 1999); DEM L50 may simply lie in 
a region of the ISM with a lower than average a/Fe ratio. 

Of course, the a/Fe ratio we observe could also result 
from inadequate X-ray models or spectra, rather than 
genuine properties of the emitting material. In X-ray 
spectral models, parameters such as temperature and 
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TABLE 11 

Hardness Ratios and Median Energies 



Region 


Q (IT fits) 


Q (2T fits) 


Q (NEI fits) 


Q (count rate) 


Median Energy (keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


D-hlvl Lou Center 


-0.67 


-0.72 


-0.63 




0.723 


DEM L50 Limb 


-0.60 a 


-0.73 a 


-0.63 a 




0.723 


DEM L50 Bright Limb 


-0.48 a 


-0.55 a 


-0.43 a 




0.767 


DEM L50 Center + Limb 


-0.63 a 


-0.72 a 


-0.63 a 


0.46 ± 0.02 


0.723 


DEM L50 SNR 




-0.56 a 


-0.54 a 


0.32 ± 0.02 


0.694 


DEM L152 Bubble 


-0.34 


-0.35 




0.75 ± 0.02 


0.869 


DEM L152 Bubble Limb 


-0.31 


-0.33 






0.854 


DEM L152 West 


-0.50 


-0.51 




0.67 ± 0.02 


0.825 


DEM L152 West Limb 


-0.48 


-0.49 






0.810 


DEM L152 West Blowout 


-0.58 


-0.59 






0.759 


DEM L152 South Bar 


-0.50 


-0.53 




0.63 ± 0.03 


0.825 


DEM L152 SNR 


-0.51 


-0.50 




0.80 ± 0.07 


0.942 



NOTE. — The hard energy band is E = 0.6 - 1.6 koV, and the soft band is E = 0.3 - 0.6 keV. 
a Using free abundance fits with the High S2 background. 



column density are degenerate with the metallicity, mak - 
ing the model difficult to constrain (Dahlem ct al. 2000). 
A variety of X-ray models may fit the same spectrum 
equally well. For instance, other authors have found that 
the need for unusual abundances in a spectral fit may 
disappear when a more complicated, mu lti-component 
X-ray spectral model is used instead (e.g. IWeaver et al.l 
[2000t IStrickland et all [20021. The ~0.12 keV resolution 
of Chandra ACIS spectra and model degeneracies make 
it unclear which precise combination of models should 
be used, however. Determining the origin of the unusual 
abundances in DEM L50 requires further investigation of 
both the LMC ISM and the X-ray models used, and the 
explanation of the observed abundances may have to wait 
for higher resolution spectra. Our current data suggest, 
however, that the non-equilibrium ionization scenario is 
the most plausible explanation for DEM L50's low ob- 
served a/Fe ratio. 

4.3.2. Luminosities 

Luminosities for all regions and models are shown in 
Table [8j with the non-equilibrium ionization fits shown 
in bold. The Limb contributes roughly 50% of the total 
emission. The Bright Limb subregion provides 13-41% of 
the total. Hardness ratios and median energies appear 
in Table [Til and indicate that the Bright Limb is clearly 
spectrally harder than the Limb as a whole, while the 
Limb and Center have comparable hardness ratios. This 
relative spectral hardness suggests that the Bright Limb 
was recently heated. For the two-temperature equilib- 
rium fit to the Limb, the Limb appears to have a softer 
hardness ratio than the Center. The two-temperature 
fit has a slightly higher absorbing column density and 
an additional lower temperature component, which ac- 
count for the difference. The temperatures of the non- 
equilibrium fits are less precise, and the larger error bars 
make the temperature differences between the regions 
unclear. The relative hardness of the Bright Limb is still 
apparent when the non-equilibrium fits are considered, 
however, and is confirmed by the median energies of the 
regions. 

As was noted in § 13.21 absorption due to intervening 
Hi may be responsible for the apparent limb-brightening 
observed in DEM L50. When the areas of the extracted 
regions are considered, the Center does appear to have 



a ~20-45% higher surface brightness than the Limb as a 
whole (Tabled]). The Bright Limb subregion, however, 
has a surface brightness comparable to or greater than 
the Center. The Hi column density has a significant ef- 
fect on the calculated luminosity, however, and as can be 
seen in Table [8j the luminosity varies substantially de- 
pending on the model. The Nn error bars in our models 
correspond to typical uncertainties in the X-ray luminos- 
ity of about 25%, but the uncertainty is up to 50% for the 
Bright Limb. Whe n Galactic Ng is i ncluded, the Njj con- 
tours observed by lOev et all (|2002l ) , shown in Figure [3J 
do agree with the modeled Nn, given the error bars in the 
fits, JVh variations within the regions, and uncertainties 
in the Nh observations. Nevertheless, we cannot conclu- 
sively determine whether the apparent brightness of the 
Bright Limb is genuine or due to line-of-sight absorption. 

The temperature and hardness of the Bright Limb sub- 
region support the idea that an off-center SNR heated 
the shell wall in that region, while the Center's bright- 
ness is most likely a result of mass-loading due to con- 
ductive evaporation from the shell walls. Mass-loading 
from swept- up clouds is unlikely given DEM L50's loca- 
tion in a void in the ISM and the lack of observed clouds 
in X-ray or other wavelengths (see Fig. [4j. We discuss 
this further in § 16.21 Metallicity enhancement is also not 
observed; while Fe is enhanced, the a elements, which 
should have the greatest effect on the cooling function, 
are not. Given the contribution of the Limb to the total 
luminosity and the lack of observed clumps in the ISM 
near DEM L50 (see Fig. gj), an off-center SNR and ther- 
mal conduction from the shell walls into the center are a 
more viable explanation for DEM L50's brightness than 
ablation or evaporation of ISM cloudlets. 

4.4. DEM L152 

Like DEM L50, DEM L152 is an X-ray-bright super- 
bubble, although as seen in Figures [5] and [3J its Hi en- 
vironment and X-ray and optical morphology differ sub- 
stantially from that of DEM L50. Consequently, DEM 
L152's derived X-ray properties also generally differ from 
DEM L50; the characteristics the two objects have in 
common may help resolve the question of why some su- 
perbubbles become X-ray-bright. We modeled the emis- 
sion from DEM L152 with one- and two-temperature 
thermal plasmas. Table [6] reports the resulting fits; the 
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column headings are the same as in Table HI except that 
the background fit, which is the same for all fits, is not 
listed. In the case of the West Limb and South Bar, it is 
not clear which fit is statistically better, so neither fit is 
shown in bold. Figure [9] shows the best fit to the Bubble 
region as an example. 

DEM L152 is a few 10 6 K hotter than DEM L50 on av- 
erage, a result of the h igher number of O stars and super - 
novae in DEM L152 (IQev fc Massevl fl995t iQevl Il996ah . 
The Bubble Limb appears slightly cooler than the Bubble 
region as a whole, which is reasonable since it makes up 
the outer edge of the Bubble and likely consists of denser, 
cooler material. The West region shows a 0.25 keV com- 
ponent as well as an extremely hot 0.79 keV compo- 
nent. This hot component appears in the West Limb 
region two-temperature fit, but not the West Blowout, 
and may also account for a substantial portion (~25%) 
of the emission from the South Bar. The localization 
of this hot component to the adjacent West Limb and 
South Bar regions suggests recent heating in this area, 
perhaps from a supernova explosion. An expanding SNR 
would also explain the limb brightening seen in the West 
region. Interestingly, iChu fe Mac Low! (|1990l ) also sug- 
gest the presence of an off-center supernova remnant in 
the sout hwest due to th e X-ra y brightness in the South 
Bar, and [Magnier et al.l (|1996l) find high-velocity gas in 
the region. 

As expected, the two blowout regions are slightly cooler 
than the Bubble region. Most of the South Bar's emis- 
sion is fit by a 0.29 keV component, while the hot compo- 
nent mentioned above accounts for the remainder. Like- 
wise, a 0.25 keV plasma characterizes the West Blowout. 
The blowouts are cooler than the 0.37 keV Bubble region 
and cooler than (or at least comparable to) the 0.33 keV 
Bubble Limb region, presumably having cooled due to 
adiabatic expansion. 

Abundances for each region show enrichment from 
core-collapse supernovae (see Table HJ) , which raises the 
a/Fe ratio. We discuss the implications of this enrich- 
ment for the expected luminosity in § 16.11 

Two clumps of X-ray emitting material, a few parsecs 
in diameter, appear in the South Bar region (see Fig. |6]). 
We extracted and fit spectra for these two knots (see 
Table [5]) , fixing the column density to the best fit value 
for the South Bar region. Due to low signal-to-noise, we 
only used single temperature fits and fixed abundances 
to solar values for the second knot. The resulting lumi- 
nosities should be viewed as estimates only. Neverthe- 
less, the luminosities are one to two orders of magnitude 
lower than the luminosity of the other regions and clearly 
do not have much effect on the total emission. We plot 
the median energy versus total count rate for the knots 
and regions of DEM L152 in Fig. [TO] The knots do not 
exhibit any detectable spectral difference relative to the 
rest of the superbubble. 

The luminosities for each region and subregion of DEM 
L152 are shown in TablefTUl the -/Vh error bars correspond 
to an uncertainty of ~T5% in the calculated luminosities. 
The Bubble region provides ~ 37% of the total X-ray 
luminosity (33% from the Bubble Limb, 4% from dif- 
fuse emission), the West region contributes ~ 41% (14% 
from the West Limb, 13% from the West Blowout, and 
14% from diffuse emission), and the South Bar provides 
~ 21% (19% from diffuse emission, 2% from the knots). 



The intrinsic surface brightnesses are unclear, as they 
depend on the model and column density assumed. For 
example, the best fits for the Bubble region imply some 
slight limb-brightening, while the two-temperature fits 
indicate the opposite. If the limb-brightening observed in 
Fig. |4] is not intrinsic, it must instead result from higher 
absorption towards the Bubble center. This hypothetical 
absorption would conflict with the Hi contours shown in 
Fig. [3l however, which show a lower column density to- 
wards the center of the Bubble. Thus, we expect that the 
observed limb-brightening in the Bubble region is real. 

A comparison of the hardness ratios and median en- 
ergies in Table [TT] shows that the South Bar and West 
Blowout regions are softer than the rest of the super- 
bubble, consistent with a blowo ut origin a n d wit h the 
previous X-ray obse rvations of IChu et all (|1993l ) and 
Ma gnier et al This variation in hardness can 

also be seen in Figure HJj. The Bubble and Bubble Limb 
regions are noticeably harder than the other regions of 
DEM L152. In general, we find softer h ardness ratios for 
DEM L152 than previous authors (e.g. iWang fe Helfandl 
119911 ). a result of fitting the column density instead of 
assuming a fixed conversion between count rate and lu- 
minosity. 

As with DEM L50, limb regions account for roughly 
half of the observed X-ray emission. The South Bar 
and West Limb regions alone provide a third of the total 
emission and may be associated with a recent off-center 
SNR as discussed above. In contrast, the two knots make 
up only 2% of the total X-ray emission, suggesting that 
emission from clumps cannot explain DEM L152's X-ray 
luminosity. 

4.5. SNRs 

The SNR N186 D appears at the northern edge of DEM 
L50. The fact that N186 D is detected in [Om] while 
DEM L50 is no t suggests th ey may be two physically 
distinct objects (Lasker Hl977D ; similar Hi kinematics in 
both objects, however, may indic ate that the SNR is as- 
sociated with the larger bubble ()Oev et al.1 [20021. and 
thus the exact relationship between the SNR and the 
superbubble remains unclear. 

SNR N186 D's X-ray spectra can be fit by a non- 
equilibrium model or a two-temperature equilibrium 
model. With free abundances, the equilibrium fits show 
a 0.19 keV component and a weaker 0.72 keV compo- 
nent. Th ese temperature s are typical of SNRs in the 
LMC (seelWilliams|lll)9^ and in our o wn Galaxy (e.g. 
iTemim et al.l 120091: |Hui fe Beckerl I2009D . If abundances 
are forced to solar values, however, the two temperature 
components are 0.08 keV and 0.29 keV. The dominant 
0.08 keV component is a lowe r tempera t ure th an any of 
the LMC remnants studied by[Williams (1999) and even 
the 0.29 keV component is on the lower end of her sam- 
ple. The free abundance fit is statistically better than 
this solar fit, with a p-value of 0.05, and shows the same 
low a/Fe ratio seen in the rest of DEM L50. Although we 
cannot rule out the possibility that this SNR is a Type 
la remnant, the consistency of the a/Fe ratio throughout 
DEM L50 suggests the abundances are associated with 
the ISM environment of DEM L50 or alternatively, that 
the SNR is not in ionization equilibrium. As discussed in 
§ 14.3.11 the non-equilibrium scenario is the more likely 
explanation. A non-equilibrium ionization fit gives an 
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Fig. 7. — The X-ray spectrum of the Center region of DEM L50. 
Table [4J. Data have been re-binned for clarity. 

ionizing timescale of 11,000 years, a temperature of 0.58 
keV, and an a /Fe ratio of .64, close to the typical LMC 
value of ~0.6 (|Smithlll999l) . 

The non-equilibrium ionization model gives a total X- 
ray luminosity in the E = 0.3-8.0 keV band of 3.73 xlO 35 
ergss -1 . The luminosity varies by an order of magnitude 
among the potential models, however, from 2 x 10 35 to 
1 x 10 36 ergss -1 , a result of the higher Nn in the mod- 
els with solar abundances. These possible luminosities 
ar e all within the range of - 10 34 -10 37 ergss" 1 found 
bv IWilliams et all (|l999l ) for SNRs in the LMC and by 
ILong et all (I2oToT) for M33 SNRs. 

The median energy for SNR N186 D, shown in Ta- 
ble [TU is softer than its hardness ratio would suggest. 
The Hi appears pa tchy in this regio n of DEM L50, and 
the observations bv lOey et al.1 (|2002[ ) show a substantial 
gradient in the Hi column density near the SNR. These 
Nn variations may be affecting our modeled tempera- 
tures in this region. 

DEM L152's SNR, 0523-679 (|Chu et al.lll99l . consists 
of a 0.38 keV plasma, close to the med ian temperatur e 
value of Williams' sample of LMC SNRs (|Williamslll999l ). 
The SNR shows a high a/Fe ratio characteristic of core- 
collapse SNe. The extremely high and unconstrained 
value of the SNR's a/Fe ratio is a result of the lower 
signal-to-noise in this small region of the bubble. The 
total X-ray lumin osity is 5.73 x 10 35 ergss -1 , typical of 
other SNRs (e.g. IWilliams et aII[T99l ILong et al]l20"l0t 
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with the one-temperature thermal plasma best fit and residuals (see 



iSeward et al.l l2010). The SNR does appear to differ spec- 
trally from the rest of the bubble, as seen in Fig. [101 

5. MODEL COMPARISONS 
5.1. Weaver model 

We modeled the central temperature and X-ray lu- 
minosity in the [0.3-2.0 keV] band for DE M L50 and 
DEM L152 using the simi larity solutions of iWeaver et akl 
(| 19771 ) as implemented bv lOev fe Massevl (|1995l ). Energy 
and mass input were not assumed to be constant; time- 
dependent energy input from stellar win ds and SNe was 
calcul ated using the mass loss rates from.Schacrc r et al.l 
( 1993) and the stella r population shown in Table IT21 (see 
lOev fe Massevlll995l for details). 

For DEM L50, we assumed an a mbient density of 1.4 
cm -3 as explained in lOevi (|1996bD . a metallicity of 0.4 
Zq and an age of 5 Myr, a t which point 2 SNe are ex- 
pected to have occurred if a lSalpeten (119551) initial mass 
function (IMF) is assumed (jOeyl Il996bl). We modeled 
DEM L152 using a density of 2.5 cm -3 (|Oev & Massevl 
1995), a metallicity of 0.4 Z Q and an age of 6 Myr, 
incorporating the effects of 4 SNe as suggested in 
lOev fe Massevb (1995) analysis of DEM L152's mass 
function. The resulting X-ray luminosities and central 
temperatures are shown in Table Q~3] and compared with 
the observed values. The observed X-ray luminosities are 
an order of magnitude higher than the standard model 
predicts. Varying the input ambient density or metallic- 
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Fig. 8. — The X-ray spectrum of the DEM L50 Center region, with the non-equilibrium ionization fit and residuals shown (see Table[5]l. 
Data have been re-binned for clarity. 



ity did increase the X-ray luminosity, but even increasing 
these parameters by an order of magnitude still resulted 
in an X-ray luminosity lower than the observed value. 

The radii calculated with the Weaver model likewise do 
not match the observed values; both bubbles have pre- 
dicted radii of more than 100 pc (130.6 pc for DEM L50, 
146.6 pc for DEM L152), whereas the observed radii are 
around 50 pc. This "growth-rate discrepancy" has been 
observed in bubbles with a range of sizes, from bubbles 
around individual Wolf- Rayet stars to the superbubbles 
we c onsider here (e.g. iSaken et al.1 Il992t iBrown et all 
ll995HOev fc Massevlll995HCappa et alj|2003t see review 
bv lOevll2009D . One possible explanation for the smaller 
observed radii is that excess radiative cooling has caused 
the bubble to lose more energy than predicted. To see if 
the high X-ray luminosities of our objects provide enough 
cooling to slow the bubble's g rowth, we use t he for mula 
for the bubble radius given bv lWeaver et al.l (|1977f ): 



R = 27n a l/5 Lli 5 tl /5 pc, 

-3 



(2) 



where no is the ambient density in cm , L^q is the me- 
chanical luminosity in 10 36 ergss -1 , and t$ is time in 10 6 
years. Assuming that supernovae dominate the mechan- 
ical energy input, we can estimate the mechanical lumi- 
nosity as the energy provided by SNe divided by the age 
of the superbubble. Two SNe and an age of 5 Myr (DEM 
L50) or four SNe and an age of 6 Myr (DEM L152) yield 
a mechanical luminosity of order 10 37 ergss -1 . The X- 



ray radiative losses for both bubbles, on the other hand, 
are on the order of 10 36 ergss -1 . Although our super- 
bubbles are X-ray-bright, the observed X-ray luminosity 
is therefore insufficient to account for the small radius 
observed in both bubbles. 

5.2. Hydrodynamic Model 

In contrast to the Weaver model, the one- 
dimensional hydrodynamic model of DEM L50 from 
lOev fc Garci'a-Segural (|2004D does correctly fit the radius 
and velocity of the bubble. The mechanical input is the 
same as before, but this model assumes a higher ambi- 
ent pressure of P/k = 1 x 10 5 cm -3 K, neglects thermal 
conduction at the shell walls, and includes the effect of a 
SNR impact on the shell. While the higher ambient pres- 
sure results in the correct radius, the X-ray luminosity 
predicted by this model is far too low, an order of magni- 
tude lower than the Weaver model, which was itself lower 
than the observed luminosity. This low luminosity sug- 
gests that thermal conduction does play a role in DEM 
L50 and is not strongly suppressed by magnetic fields. 
If thermal conduction is occurring, the bubble temper- 
ature should be lowest at the outer edge, where dense, 
cold evaporated material is present, and should increase 
toward the bubble center. The equilibrium fits to the 
DEM L50 Limb and Center regions do show this tem- 
perature contrast; the effective temperature of the Limb 
is only 0.09 keV, while the value increases to 0.15 keV 
for the Center region. The non-equilibrium fits, on the 
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Fig. 9.— The X-ray spectrum of the Bubble (Shell 1) region of DEM L152 
residuals (see Table [6}. Data have been re- binned for clarity. 
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with the one-temperature thermal plasma best fit and 
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Fig. 10. — Median Energy vs. Total Count Rate for regions in 
DEM L152. The plot shows rough spectral differences between 
the regions. Error bars were estimated using the quanti l e erro r 
estimation approach outlined in Appendix B of Hong et al.| 1120041 ). 

other hand, suggest that the Limb is slightly hotter than 
the Center. The possibility of a cooler Limb is within the 
error bars, however, and the contribution of the heated 
Bright Limb region may raise the average temperature 



of the Limb. 

The X-ray luminosity at different times is shown in Ta- 
ble Q31 and compared with the Weaver model luminosity 
for DEM L50 and the observed luminosity. A super- 
nova goes off at t—5.26 Myr, and the SNR increases the 
bubble's luminosity by several orders of magnitude. The 
emissivity profiles of the bubble before and after the SNR 
hits the shell walls are shown in Figure [TT] The SNR im- 
pact enhances the emissivity of the bubble's limb and also 
increases the shell velocity from 1 kms" 1 to 23 kms -1 in 
agreement with the observed value. Although the lumi- 
nosities of this model are unrealistic, the response of the 
luminosity and shell velocity to the SNR impact suggests 
that SNR impacts may be important in our superbub- 
bles. 

Neither the Weaver model nor the 1-D hydrodynamic 
model is able to match the observed luminosities of the 
superbubbles. The hydrodynamic model suggests, how- 
ever, that mass-loading via thermal evaporation from the 
shell walls must play a role in these bubbles. Mass- 
loading from the shell walls and heating from SNR im- 
pacts are consistent with our observations and could ac- 
count for most of the superbubbles' X-ray emission. 

6. ADDITIONAL X-RAY ENHANCEMENT MECHANISMS 

6.1. Metallicity 

Silich] ()200lD suggests that since emission lines of a 
elements dominate the soft X-ray emission, metallicity 
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TABLE 12 

Stellar Populations 



Bubble 
(1) 


85 M 
(2) 


60 Mq 
(3) 


40 M Q 
(4) 


25 M 
(5) 


20 Mq 
(6) 


15 Mq 

(7) 


12 Mq 

(8) 


DEM L50 
DEM L152 


1 
1 


1 

3 


3 
1 


1 

11 


8 


6 
12 


9 



Note. — Number of stars of a giv en mass assumed for the Weaver-based models of 
DEM L50 and DEM L152. Data from lOev fc Massevl 1 11999 1 and lOevl dl996bl1 . 



TABLE 13 

Weaver Model Results 



Bubble Model 


no 


Age 


L x 


[0.3-2.0 kcV] 


kT c 


[0.3-2.0 keV],obs a 


fcT c ,obs 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


DEM L50 Weaver 


1.4cm -3 


5 Myr 


1.31 


xlO 35 ergss -1 


0.27 keV 


2.0 - 4.5 X 10 36 ergss" 1 


0.15 keV b or 0.28 keV c 


DEM L152 Weaver 


2.5 cm -3 


6 Myr 


5.05 


xlO 35 ergss -1 


0.32 keV 


5.4- 5.7 x 10 36 ergss -1 


0.37 keV d 



a Excludes SNRs. 

k From single-temperature equilibrium fit to Center region. 
c From non-equilibrium fit to Center region. 
^ From best fit to Bubble region. 



TABLE 14 

Hydrodynamic Model Luminosities 



Energy Band 


t =5.152 Myr 


t =5.264 Myr 


t =5.376 Myr 


t =5.488 Myr 


t =5.600 Myr 


L x Weaver 


L x Obs a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0.3-2.0 keV 


3.7 x 10 31 


1.2 x 10 34 


8.9 x 10 32 


6.0 x 10 33 


3.0 x 10 33 


1.3 x 10 35 


2.0 - 4.5 x 10 36 


2.0-8.0 keV 


1.3 x 10 31 


2.5 x 10 33 


1.1 x 10 33 


5.3 x 10 32 


3.9 x 10 32 




2.8 x 10 34 



Note. — Results refer to DEM L50 only. 
a Excludes SNRs. Luminosities in units of ergss -1 . 

enhancement from SNe could explain the anomalously 
high X-ray luminosities observed in many superbubbles. 
As discussed earlier (§5 l4T3P~4l) . DEM LI 52 shows an a 
enhancement, but DEM L50 may not, indicating that a 
different explanation would be required for that bubble. 

Assuming the init ial stellar population proposed by 
IQev fc Massevl (|1995|) , approximately four SNe have gone 
off in DEM L152, from one 85 Mq progenitor and three 
60 Mq progenitors (see Table [T^l). Th ese SNe should in- 
put a total 65 Mq of O (lMaeder|[1992ft an d 0.2-0.8 Mq of 
Fe (see discussion in ISilich et al.l 1200 jT) . Assuming a 
number density of 0.01-0.2 cm~ 3 , spherical geometry, 
a radius of 50 pc, and 10% He abundance, the total 
mass enclosed by the bubble should be between 100 and 
3000 Mq. As the inferred number density from the ob- 
served spectrum (see Equation[T]) is 0.2 cm -3 , we assume 
an enclosed mass of 300 Mq in what follo ws. We use 
Equations 9 and 10 from ISilich et al.l (|2001[) to calculate 
the resulting metallicity: 



Zo = 



M e j y o/Z so \,o + ZismM cv 



M, 



Me 



Zfc 



-^ej,Fe/%sol,Fe + ZismM c 



Mr. 



Me 



(3) 



(4) 



M c j is the mass ejected in SNe, M ev is the mass evapo- 
rated from the shell walls, and Z S oi,o and Z S oi,Fc are the 



solar abundances by mass from Greves se et al.l (fl 9961. 
Z so1 .q = 0.0083 and Z so i iPe = 0.00126. iWeaver et al.l 
(|1977f) predict that the evaporated mass will be the dom- 
inant source of mass in the bubble interior. We assume 
that, apart from the mass ejected by the four SNe, the 
remainder of the mass was contributed by evap oration. 
Assu ming stellar remnant masses of ~ 2 Mc?> (iMaederl 
119921) and the initial masses from lOev k, Massevl (|1995f)~ 
the total mass returned to the ISM by the four most mas- 
sive stars is 257 Mq, resulting in a M ev of ~ 2700 Mq. 
Then, from equations [3] and |4j the predicted metallicities 
are Z = 3 Z Q and Z Fc = 0.4 - 0.6 Z Q . If the bubble 
contains less evaporated mass, these metallicities will be 
higher. 

As noted by ISilich et al.l ((2001), the X-ray emissivity 
is proportional to the metallicity. Assuming the a el- 
ements dominate the soft X-ray emission, the factor of 
10 increase from Zlmc of 0.4 Zq to Zq of 3 Zq would 
be sufficient to explain the factor of 10 increase in DEM 
L152's luminosity (see Table ITBI) This predicted metal- 
licity does not match the observed a abundance in DEM 
L152, however (see Table[6]), which is subsolar in all cases 
with reasonable error bars. The a/Fe ratio is a more 
reliable diagnostic of the metallicity, and it is possible 
that the a abundance from the spectral fits is incorrect. 
Given the additional uncertainties in the bubble's den- 
sity, geometry, and resulting mass and in the number of 
SNe that have occurred, this discrepancy between the 
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Fig . 11. — Emission measure along the line of sight for the 1-D hy - 
drodynamic models of DEM L50 bv lOev k, Garci'a-Segural H2004h . 
The top panel shows the modeled bubble at t=5.152 Myr, before 
the SN at 5.26 Myr goes off. The subsequent panels show the 
bubble after the SNR expands into the shell walls. 

predicted and observed metallicities is not unexpected. 
The predicted ct/Fe ratio from the 4 SNe is 5-8, while 
Table |6] shows typical a/Fe ratios of around 3. When 
errors are taken into account, the observed a/Fe ratio 
could be consistent with the predicted value, and there- 
fore it is possible that the higher a abundance is also 
present. 

If we use the observed a abundances, which have an 
average value of ~ 0.6 Zq and are generally consistent 
with 1 Zq , the enhancement over the Weaver model only 
accounts for 14-20% of the observed emission. At most, 
the observed a abundance can account for the diffuse 
emission in the Bubble and West regions and the West 
Blowout, but not for the emission from the Bubble Limb, 
West Limb, or South Bar. If, on the other hand, we ac- 
cept that this observed a abundance is inaccurate and 
use the predicted Zo, the enhanced metallicity explains 
the X-ray luminosity but not the limb-brightened fea- 
tures or enhanced expansion velocities. 

A similar analysis for DEM L50 yields a predicted Zo 
of 5Zq and predicted a/Fe of 7-13, where we have as- 
sumed M eji0 = 36 M Q , M cj!Fo = 0.1 - 0.4 M Q , a bubble 
radius of 50 pc, and njj = 0.06 cm -3 , derived from the 
non-equilibrium fit normalization. The predicted metal- 
licities are much higher than the approximately solar val- 
ues obtained from the non-equilibrium fit, and as with 



DEM LI 52, the observed metallicities cannot account for 
the majority of the X-ray enhancement. 

6.2. Mass-loading by clouds 

Another possible explanation for the high X-ray lumi- 
nosities of DEM L50 and DEM L152 is mass loading from 
the evaporation or ablation of overrun clouds. While the 
distribution of mass from evaporation at the shell walls 
should be uniform over relatively large scales, evaporat- 
ing clouds should appear as small clumps of emission in 
our data. DEM L50 is located in a void and does not ex- 
hibit any noticeable dumpiness in its emission. The two 
knots we observed in DEM L152 and its more crowded HI 
environment, on the other hand, show that mass-loading 
by clu mps could be occurring in that bubble. iSilich et "all 
(1996) investig ate clouds wit h n e ~ 10 cm -3 and radii of 
3-5 pc, while lOrlando et ail (|2005l) find their modeled 
shocked clouds have number densities of 4 cm -3 , radii 
of 1 pc, and are surrounded by a corona with density 
0.4 cm -3 . To encompass this range of possible cloud pa- 
rameters, we calculated the emission from clouds with 
each of these three densities (n e =0.4, 4, and 10 cm -3 ), 
a temperature of 10 6 K, a metallicity of 0.4 Zq, and radii 
of 0.5 pc, 1 pc, and 5 pc, all of which should be resolved at 
Chandra's ~ 1" resolution. The X-ray emissivity in the 
0.2-4.5 ke V band is approxim ately A x — 9 x 10 -24 x Z erg 
s -1 cm -3 (Sili ch et aflll996D ; this emissivity will result 
in more emission than we expect to see in the 0.3-2.0 keV 
band, but it should allow us to roughly estimate the type 
of clouds that could be det ectable. Assuming a median 
photon energy of 0.55 keV (|Orlando et al.ll2010D . we can 
calculate the predicted count rate per arcsec 2 (Table [15]). 
A comparison with the X-ray contours in Figures [5] and [3] 
shows that the knots with radii of 0.5 pc are potentially 
detectable if they have a density of ~ 10 cm -3 , while no 
knots with a density of ~ 0.4 cm -3 are detectable. Inter- 
estingly, the detected "knots" in DEM L152 have a count 
rate of approximately 6.4 x 10 -6 s -1 arcsec -2 , which is 
within the range of values in Table [T5] and a density of 
~ 1 cm -3 , found using Equation [TJ assuming a radius of 
2.5 pc. These values suggest that these knots may indeed 
be genuine clouds swept up by the bubble. 

Another possibility is that some clumps appeared as 
point sources in the images (Fig. [T]). We would expect 
predominantly soft emission from overrun clumps, while 
other point sources, such as background AGN or X-ray 
binaries, should appear in the hard band as well. Only 
two point sources, one in DEM L50 and one in DEM 
L152, appear in the soft band alone. The point source 
in DEM L50 is quite faint, containing only ~ 15 counts, 
and its identity is unclear. The soft point source in DEM 
L152 is located outside the main bubble. This source is 
likely associated with the highly ionized N44 C region 
rather than an independent cloud swept up by DEM 
L152. In any case, the lack of multiple soft point sources 
suggests that dense, overrun globules are not present to 
any significant degree in either bubble. 

Although DEM L50 and DEM L152 may have over- 
run additional clouds which are below our detection 
limit or already destroyed, their morphology and spec- 
tral properties suggest that mass-loading by clouds is 
not dominating the bubbles' evolution. Modeled bub- 
bles with high levels of mass loading tend to show a 
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TABLE 15 
Predicted Cloud Count Rates 



ra e (cm 3 ) 


i?, c l=0.5 pc 


i? cl =1.0 pc 


i?, cl =5.0 


pc 


(1) 


(2) 


(3) 


(4) 




0.4 


2.5 x 10~ 9 


5.0 x 10~ 9 


2.5 x 10 


-8 


4.0 


2.5 x 10" 7 


5.0 x 10~ 7 


2.5 x 10 


-6 


10.0 


1.5 x 10" 6 


3.1 x 10~ 6 


1.5 x 10 


-5 


Note. — 


Count rates 


are given for 


the 0.2-4.5 


keV 


energy band 


in photons s — 1 


arcsec - 2 . 







smoother distribution of density and temperature across 
the bubble, as the clouds add cooler gas to the sys - 
tem (see e.g. lArthur fc Hennevl [l99l IPittardl f2007h 
The increased central density and reduced central tem- 
perature may result in a flat surface brightness profile 
(lArthur fc Hennevl 1996) or even a centrally-peaked pro- 
file (jSilich et al.l 119961) . A limb-brightened morphology, 
on the other hand, appears correlated with supernova 
activity. No pre-supernova superbubbles yet studied 
show limb-brightening, while this mo rphology is com - 
mon in post-supernova superbubbles (|Chu et al.l 120031 ) . 
DEM L50 and DEM L152's limb-brightened morpholo- 
gies likely reflect their supernova history rather than a 
history of high mass-loading. 

Mass-loading sho uld reduce the cen tral temperature 
of a bubble (e.g. iSilich et al.l I1996Q . whereas DEM 
L152's central temperature is quite high (~ 0.37 keV), 
higher than predicted by the Weaver model. DEM L50's 
lower central temperature of 0.15 keV in the equilib- 
rium fit could be consistent with mass-loading, although 
uncertainties in the number of SNe or energy input 
in DEM L50 could also account for the slightly lower 
temperature. The temperature derived from the non- 
equilibrium fit to DEM L50 is completely consistent with 
the Weaver model predictions, however, requiring no ad- 
ditional mass-loading. Mass-loading by overrun clouds 
or evaporation from the shell walls is likely occurring to 
some extent in both bubbles, but it does not fully explain 
the bubbles' X-ray emission or large-scale properties. 

Dense cloudlets have been detected in other superbub- 
bles, but their presence does not necessarily imply a high 
X-ray luminosity. For instance, although clouds are de- 
tected in the superbubble DEM L192 and the bubble ap- 
pears to be expanding into an inhomogeneous ISM, its X- 
ray luminosity in the 0.3-3.0 keV band is o nly moderate, 
~ 4 x 10 35 ergss" 1 (|Cooper et al.ll200l . iCooper et all 
( 2004) note a correlation between the Ha and X-ray sur- 
face brightnesses in this bubble, however, and suggest 
that mass-loading is occurring through evaporation at 
the shell walls. This analysis is consistent with our find- 
ings in this section and in § 15.21 thermal conduction and 
the subsequent evaporation of shell material is important 
in DEM L50 and DEM L152, while the contribution of 
swept-up clouds to the X-ray emission is negligible. 

Thus, we see evidence that both metallicity enhance- 
ment and mass-loading by clouds has taken place in these 
superbubbles. The metallicity enrichment is too low to 
account for the observed luminosities, however. Likewise, 
few candidate clouds are observed in the superbubbles 
and the X-ray emission is inconsistent with a high mass- 
loading scenario. Metallicity and mass-loading by clouds 
do not completely explain the observed X-ray luminosi- 



ties. 

7. CONCLUSIONS 

We present high spatial resolution Chandra ACIS-S im- 
ages and spectra of the diffuse emission in two X-ray- 
bright LMC superbubbles, DEM L50 and DEM L152. 
These objects have X-ray luminosities on the order of 
10 36 ergss -1 , an orde r of magnitud e high er than pre- 
dicted by the standard I Weaver et al.l (|1977|) wind-blown 
bubble model. The level of detail in the Chandra ob- 
servations allows us to evaluate the relative importance 
of the various mechanisms that contribute to the X-ray 
luminosity by constraining the spatial origin of the emis- 
sion and revealing X-ray spectral variations across the 
superbubbles. 

With these high-quality images, we isolated the diffuse 
X-ray emission by excluding point sources and examined 
the spectra of separate regions within each superbubble. 
X-ray images of both bubbles show limb-brightened mor- 
phologies. Most of the limb-brightening in DEM L50 is 
the result of a higher central absorption column density, 
while the limb-brightening at the bubble's southern end 
may be real. Likewise, the intervening absorption does 
not account for the limb-brightening in Shell 1 of DEM 
L152. Hardness ratios indicate that these same regions, 
the Bright Limb subregion of DEM L50 and Shell 1 of 
DEM L152, are spectrally harder than the other regions 
in the superbubbles. 

Spectral fits to the superbubbles account for the varia- 
tions in absorption column density from region to region 
and reveal the temperature differences and luminosity 
contribution of each region. The spectra of both objects 
show evidence of SNR impacts on the shell walls. In par- 
ticular, DEM L50's Bright Limb has a harder spectrum, 
higher temperature, and higher surface brightness than 
the rest of the bubble. These spectral properties support 
the idea that a SNR recently heated the Limb region. 
The limb-brightened Shell 1 region in DEM L152 also has 
a higher temperature and surface brightness compared 
to other regions of the superbubble. A high temperature 
component in the limb-brightened West region and the 
adjacent South Bar blowout may be another indicator of 
recent heating due to SNe. 

Our spectral fits also revealed unexpected abundance 
ratios in DEM L50, which suggest the superbubble is 
not in collisional ionization equilibrium. DEM L50's low 
a/Fe ratio is the opposite of the expected a-enrichment 
from core-collapse SNe. While we cannot rule out the 
possibility of a Type la SN in the area, the consistency 
of the ratio across the bubble suggests that either the 
bubble happens to lie in an unusual ISM environment 
or, more likely, DEM L50 is not in ionization equilib- 
rium. Assuming a non-equilibrium situation, the derived 
ionization timescales are on the order of 10 5 years, which 
is consistent with the other indications of a recent super- 
nova in the bubble. 

We compared our observations with both the stan- 
dard wind-blown bubble model (We aver et al.l I1977D 
and a one-dimensional hyd rodynamic model based on 
IQev fc Gar cia-Scgura (2004). Both models include the 
time-dependent energy input of the stellar population. 
As shown by the Weaver model predicts a total X-ray 
luminosity an order of magnitude below the observed 
value. While the 1-D hydrodynamic simulation repro- 
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duces the size and expansion velocity of DEM L50 by 
assuming a higher ambient pressure and including the ef- 
fects of SNR impacts, the neglect of thermal conduction 
produces X-ray luminosities far below our observed lu- 
minosities. This discrepancy indicates that thermal con- 
duction is important in these bubbles, as evaporation 
of material from the shell walls enhances the density of 
the interior. The SNR impact does significantly brighten 
the modeled bubble, particularly at the limb, which sup- 
ports the off-center SNR ex planation for the existe nce of 
X-ray-bright superbubbles (|Chu fc Mac Lowlll99Cl[ ). 

Our observations do not reveal signs of significant 
metal enrichment, an a lternative explana tion for the high 
observed luminosities (Sili ch et al.ll200il ). A higher pro- 
portion of metals could enhance the X-ray emissivity, but 
neither superbubble exhibits a sufficiently high metal- 
licity to account for the anomalous X-ray luminosity. 
Abundances determined from X-ray observations are un- 
certain, however, and the bubbles may be more enriched 
than the spectral fits suggest. While a high metallic- 
ity could raise the luminosity, it would still not explain 
the observed morphologies or large expansion velocities 
of the objects. 

We consider the hypothesis that mass-loading through 
the ablation or evaporation of swept-up clouds could in- 
crease the density of the superbubbles and raise the X- 
ray l uminosity (e.g. iMcKee et al.l 119841 : iPittard et al.l 
|2001[ ). With Chandra's ~ 1" resolution, we expect to 
detect clouds a few parsecs in diameter. Few candidate 
clouds are observed, and they contribute negligible emis- 
sion to the superbubbles. Although some mass-loading 
is likely present and the superbubbles may have overrun 
more clouds in the past, the effect of clouds on the X- 
ray emission is insignificant compared to the luminosity 
enhancement from SNR impacts. 

An examination of the energy budget for both super- 
bubbles further clarifies the role of off-center SNRs. The 
relatively small Bright Limb subregion of DEM L50 gen- 
erates roughly a quarter of the total luminosity, while 
the limb as a whole accounts for half of the total X-ray 
luminosity. In DEM L152, the Shell 1 Limb and West 
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Limb regions alone provide about half of the total X-ray 
luminosity. All these regions exhibit signs of supernova 
activity. On the other hand, the contribution of metal- 
licity enrichment and mass-loading from clouds to the 
X-ray luminosity is much less. Metallicity enrichment 
accounts for at most 20% of the excess emission, while 
dense clouds account for less than 2% of the total emis- 
sion. SNR impacts and evaporation of material from the 
shell walls dominate the observed X-ray luminosities. 

Our Chandra images and spectra of DEM L50 and 
DEM L152 are the first observations of X-ray-bright su- 
perbubbles at a sufficiently high spatial and spectral res- 
olution to distinguish between various emission mecha- 
nisms. While mass-loading and metallicity enrichment 
due to SNe must be occurring to some degree, our ob- 
servations confirm that thermal evaporation from the 
shell walls and off-center SNR impacts are the primary 
sources of the enhanced X-ray luminosity. SNR impacts 
increase the shell expansion velocity, lead to higher radia- 
tive losses, and may contribute to blowouts and the sub- 
sequent mixing of ISM and hot, enriched gas. Although 
the X-ray luminosities are not high enough to explain 
the superbubble growth-rate discrepancy, X-ray radia- 
tive losses will decrease the energy available for transport 
to the ISM. Unlike the cloud mass-loading scenario, SNR 
impacts and thermal conduction should affect most su- 
perbubbles at some point in their evolution, irrespective 
of their environment. Off-center SNR impacts in super- 
bubbles play a critical role in the evolution of the ISM by 
enhancing the mixing of hot, enriched gas with ambient 
material, while at the same time leading to X-ray cooling 
that decreases the available energy. 
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